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Abstract 

- - The  compact  range  allows  for  the  measurement  of  radar  cross  section  and 

antenna  patterns  in  a  relatively  small  room.  An  offset-feed  parabolic  reflector  is 
normally  used  to  approximate  a  plane  wave  in  the  target  zone;  however,  some  of  the 
energy  transmitted  by  the  feed  will  strike  the  ceiling,  walls,  and  floor  of  the  room. 
These  stray  signals  then  scatter  in  all  directions.  To  reduce  the  level  of  scattered 
energy,  the  room  is  typically  lined  with  pyramidal-shaped  radar  absorbing  material. 
However,  the  behavior  of  this  material  is  not  well  understood. 

A  diffraction  formulation,  based  on  the  Uniform  Geometrical  Theory  of  Diffrac¬ 
tion  is  developed  for  a  lossy  dielectric  corner.  Using  this  formulation,  a 
computer  code  is  written  that  calculates  the  bistatic  scattering  from  a  pyrami¬ 
dal  absorber  tip.  Sample  results  display  some  features  of  scattering  from  a  single 
pyramid.  Calculations  are  then  comj>ared  with  backscatter  measurements  of  a 
single  pyramid,  and  with  bistatic  measurements  fSf'of  an  ubsorber  wall.  Next,  a 
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general  purpose  computer  code  is  written  which  calculates  the  scattering  into  the 
target  zone  of  a  compact  range  from  the  pyramidal  absorber  which  lines  the  room. 
(A  focus-fed  semi-circular  parabolic  reflector  with  rolled  e  dge  and  a  skirt  '(4]’illu- 
minates  the  room.  ample  calculations  show  the  total  power  scattered  into  the 
target  zone,  and  how  strongly  various  locations  in  the  room  contribute  to  this  total 
scattered  power.  Finally,  calculations  ,are  compared  with  experimental  data  [5]  in 

,  y*  (rk* j ,) 

which  a  6  x  6  aluminum  plate  (mounted  in  the  Ohio  State  University  compact 
range  target  zone)  was  used  to  measure  bistatic  scattering  from  the  pyramidal 
absorber  on  the  ceiling. 
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CHAPTER  I 


INTRODUCTION 

>  1.1  Statement  of  the  Problem 

The  measurement  of  radar  cross  section  (RCS)  or  of  an  antenna  pattern  re¬ 
quires  that  the  target  or  antenna  being  tested  be  illuminated  by  a  uniform  plane 
•wave.  This  plane  wave  is  approximated  in  a  far-field  test  range  simply  by  placing 
the  test  object  in  the  far  zone  of  the  transmitter,  where  the  spherical  wavefront  of 
the  transmitted  field  can  be  considered  approximately  planar.  Depending  on  the 
size  of  the  test  object  and  the  wavelength  being  used,  this  far  zone  requirement 
typically  leads  to  an  outdoor  facility.  The  compact  range  offers  an  alternative 
to  the  far-field  test  range  for  the  measurement  of  RCS  or  antenna  patterns.  A 
compact  range  normally  uses  an  offset-feed  parabolic  reflector  to  approximate  a 
plane  wave  in  a  relatively  small  chamber.  The  small  size  of  a  compact  range  makes 
it  practical  as  an  indoor  facility,  thus  free  from  adverse  weather,  interference,  or 
external  monitoring.  Of  course,  not  all  of  the  energy  transmitted  by  the  feed  will 
i  be  directed  into  the  desired  plane  wave.  Rather,  some  will  strike  the  ceiling,  walls, 

and  floor  of  the  room.  These  stray  signals  then  scatter  in  all  directions.  They 
degrade  the  nature  of  the  incident  plane  wave  in  the  target  region,  and  they  also 
return  directly  to  the  receiver,  corrupting  the  desired  return  from  the  target.  Thus, 
they  affect  the  accuracy  of  the  range  and  the  size  of  its  “sweet  spot"  (the  region 
where  the  illumination  satisfactorily  approximates  a  plane  wave).  One  method  to 
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eliminate  these  stray  signals  might  be  to  make  the  room  large  enough  so  that  they 
could  all  be  time-gated  out;  however,  this  would  be  prohibitively  expensive.  In 
a  room  of  practical  size,  some  of  these  unwanted  signals  can  be  time-gated  out, 
others  cannot.  To  reduce  the  level  of  scattered  energy,  the  room  is  typically  lined 
*  with  pyramidal-shaped  radar  absorbing  material.  However,  the  behavior  of  this 

pyramidal  absorber  is  not  well  understood,  especially  for  compact  range  applica- 
^  tions. 

In  this  study,  a  scattering  formulation  based  on  the  uniform  geometrical  theory 
of  diffraction  (UTD)  is  developed  to  describe  the  fields  diffracted  by  a  pyramidal 
absorber  tip.  This  is  then  used  to  generate  a  computer  code  for  the  calculation 
of  bistatic  scattering  from  the  tip  of  an  absorber  pyramid.  With  this  code,  and  a 
reflector  analysis  code  that  computes  the  fields  incident  on  the  pyramidal  absorber 
in  the  chamber,  a  general  purpose  computer  code  is  developed  which  calculates  the 
energy  scattered  by  the  pyramidal  absorber  to  the  target  region  of  the  compact 
range  chamber.  Finally,  various  scattering  measurements  are  made  to  verify  the 
calculated  results.  These  range  from  the  backscattcr  of  a  single  pyramid  to  the 
absorber  scatter  directed  into  the  compact  range  target  zone  from  the  ceiling. 

1.2  Motivation 

Pyramidal  absorber  material  lias  been  widely  used  for  many  years  to  reduce 
't  the  clutter  level  in  anechoic  chambers.  However,  as  stated  above,  the  behavior 

of  pyramidal  absorber  material  is  not  well  understood.  The  actual  levels  of  stray 
signals  scattered  from  the  absorber  in  a  compact  range  chamber  have  not  previously 
been  calculated.  Such  information  would  prove  useful  in  the  design  of  a  compact 
range  chamber,  and  thus  motivates  this  research  effort. 


1.3  Approach 


The  first  goal  in  this  study  is  the  development  of  a  scattering  formulation  for 
a  (lossy)  dielectric  pyramidal  tip.  The  pyramidal  tip  is  analyzed  as  a  corner  in  a 
three-dimensional  wedge,  which  has  four  planar  faces  and  four  edges  intersecting 
at  a  common  point.  Thus,  a  dielectric  corner  diffraction  formulation  is  needed. 

In  Chapter  II,  the  expressions  for  diffraction  from  perfect  conductors  are  re¬ 
viewed.  These  include  two  and  three-dimensional  wedge  diffraction  [1,2, 3, 4],  along 
with  corner  diffraction  [5].  An  understanding  of  these  is  necessary  for  one  to  ex¬ 
tend  the  perfectly-conducting  corner  diffraction  solution  to  the  case  of  a  dielectric 
corner. 

Chapter  III  begins  by  investigating  dielectric  modifications  that  have  been 
made  to  the  UTD  for  simpler  geometries.  The  first  was  developed  by  Burnside  and 
Burgener  [6],  who  studied  the  problem  of  edge  diffraction  from  a  two-dimensional 
dielectric  slab.  This  was  based  on  the  effect  of  the  dielectric  on  the  discontinuities 
at  the  shadow  boundaries,  and  accounted  for  the  multiple  interactions  occurring 
within  the  dielectric  slab.  Next,  DeWitt  [7]  considered  a  similar  modification 
for  the  more  general  problem  of  a  two-dimensional  dielectric  wedge.  Owing  to  the 
added  complexity  of  the  wedge  problem,  the  method  of  solution  was  found  through 
a  heuristic  argument.  However,  comparisons  with  a  Neumann  series  solution  by 
Rawlins  [8]  show  good  agreement. 

Chapter  III  then  proceeds  to  further  thes<  efforts.  Sections  3.5  and  3.6  extend 
the  dielectric  modification  of  [7]  to  the  problems  of  three-dimensional  dielectric 
wedge  and  corner  diffraction,  respectively,  The  most  significant  complication  here 
involves  the  matter  of  field  polarization  inherent  in  the  three-dimensional  problem. 
An  interesting  result  is  that  the  dielectric  material  causes  a  cross-coupling  in  the 


polarizations  of  the  incident  and  diffracted  fields. 

Chapter  IV  applies  the  UTD  dielectric  corner  diffraction  solution  to  the  prob¬ 
lem  of  the  absorber  pyramid.  The  interpretation  and  calculation  of  the  variables  in 
the  dielectric  corner  diffraction  solution  represents  a  significant  task,  and  is  thus 
described  in  detail.  The  resulting  code,  referred  to  as  the  Absorber  code,  com¬ 
putes  the  bistatic  scattering  from  a  pyramidal  absorber  tip.  Sample  calculations 
for  various  monostatic  and  bistatic  geometries  are  then  presented. 

Chapter  V  considers  the  absorber  scattering  that  takes  place  in  a  compact 
range  chamber.  The  fields  from  the  range  reflector  which  illuminate  the  absorber 
are  calculated  through  a  modified  version  of  the  Semi-Circular  Compact  Range 
Reflector  Code  [9],  while  the  fields  subsequently  scattered  by  the  absorber  to  the 
target  region  are  calculated  by  means  of  the  Absorber  code.  (This  is  done  in  a 
general  purpose  computer  code,  program  CHAMBER,  which  calls  the  reflector  and 
absorber  codes  as  subroutines.  CHAMBER  allows  many  parameters  to  be  varied, 
and  provides  several  output  options.)  The  method  of  calculation  is  first  described 
in  some  detail.  Sample  calculations  of  absorber  scattering  in  a  compact  range  arc 
then  presented. 

Chapter  VI  provides  experimental  verification  of  absorber  scatter  calculations. 
First,  the  backscatter  of  a  single  pyramidal  tip  is  measured,  and  is  seen  to  compare 
well  with  calculations.  The  topic  of  addition  of  the  scattered  fields  from  several 
pyramids  is  then  discussed,  since  this  is  relevant  to  the  remaining  experiments. 
Next,  calculations  are  compared  with  measurements  of  bistatic  scattering  from  a 
wall  of  pyramidal  absorber.  The  results  are  very  encouraging,  since  the  calculations 
predict  both  the  magnitude  and  the  frequency  behavior  of  the  measured  absorber 
scatter.  Finally,  an  experiment  is  described  in  which  the  reflected  fields  from  a 
large  plate  are  directed  toward  the  ceiling  absorber  (of  the  Ohio  State  University 


compact  range  chamber)  and  are  then  bistatically  scattered  by  the  absorber  to  the 
reflector  and  received  by  the  compact  range.  Again,  calculations  agree  well  with 
measurements. 

Chapter  VII  presents  a  summary  and  conclusions.  In  appendix  A,  the  opera¬ 
tion  of  program  CHAMBER  is  discussed. 


CHAPTER  II 


DIFFRACTION  FROM  PERFECT  CONDUCTORS 

2.1  Introduction 

This  chapter  begins  by  describing  the  UTD  solution  for  scattering  from  a 
perfectly-conducting  two-dimensional  wedge.  Two  ray-fixed  coordinate  systems 
are  then  introduced  and  used  to  put  the  UTD  solution  for  a  perfectly-conducting 
three-dimensional  wedge  into  a  compact  form.  Lastly,  diffraction  coefficients  for  a 
corner  in  a  conducting  flat  plate  and  a  three-dimensional  wedge  are  presented. 

2.2  Two-Dimensional  Wedge  Diffraction 

Consider  an  electric  or  magnetic  line  source  illuminating  an  infinite,  perfectly- 
conducting  wedge  formed  by  two  plane  surfaces,  as  shown  in  Figure  1.  The  total 
UTD  solution  for  this  problem  consists  of  an  incident,  reflected,  and  diffracted 
field.  The  incident  field  is  the  source  field  in  the  absence  of  the  wedge  (although 
one  must  take  into  account  the  shadowing  of  the  wedge).  The  reflected  field  is 
that  field  reflected  from  the  surface  of  the  wedge  with  the  edge  ignored.  Region  1 
contains  incident,  reflected,  and  diffracted  rays;  region  II  contains  incident  and 
diffracted  rays;  and  region  III  contains  only  diffracted  rays.  The  incident  and 
reflected  fields  constitute  the  classical  geometrical  optics  (GO)  solution.  Alone, 
this  solution  predicts  zero  fields  in  region  Ill  (the  shadow  region),  and  disconti¬ 
nuities  exist  at  both  the  incident  and  reflected  field  shadow  boundaries.  It  is  the 


Field 


Figure  1:  UTD  ray  paths  of  a  line  source  radiating  in  the  presence  of  a  wedge. 

diffracted  held  that  exi6t6  in  the  shadow  region,  and  it  is  the  diffracted  held  that 
combines  with  the  geometrical  optics  held  to  produce  a  total  held  which  is  contin¬ 
uous  across  both  shadow  boundaries.  Keller  [1,2,3}  hrst  introduced  the  diffracted 
field  in  his  geometrical  theory  of  diffraction  (GTD),  yielding  a  simple  geometric 
interpretation  of  the  fields  associated  with  the  currents  excited  at  the  edge  of  the 
wedge.  However,  one  limitation  of  his  solution  was  that  it  was  singular  at  the 
shadow  boundaries,  where  the  diffracted  held  is  most  significant.  Kouyoumjian 
and  Pathak  [4],  using  asymptotic  methods  and  the  exact  eigenfunction  solution 
for  the  wedge,  developed  a  uniform  geometrical  theory  of  diffraction  aiisatz  within 
the  context  of  Keller's  solution.  Their  solution  eliminates  Keller's  singularity  by 
introducing  a  multiplicative  factor  that  approaches  zero  at  the  shadow  boundaries. 
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The  total  UTD  solution  for  the  wedge  is  given  by 


utotal  -  uinc  +  uref  +  udif  (2.1) 

where  u  represents  an  electric  (magnetic)  scalar  field  for  illumination  by  an  electric 
(magnetic)  line  source.  The  line  source  is  located  at  (p' ,<t>')  with  respect  to  the 
edge,  while  the  receiver  is  located  at  (p,</>).  The  incident  field  is  given  by 

I  Ke  yr*  in  regions  I  and  II,  and 

u  "“  =  (2.21 


y  0  in  region  III 

where  p{  is  the  distance  from  the  source  to  the  receiver,  Jfe  is  the  wavenumber  of 
the  medium  surrounding  the  wedge,  and  K  is  a  complex  constant.  The  reflected 
field  is  given  by 

„«/  =  {  iKt:jzr  *e>»»  ,,  3, 


(  0  in  regions  II  and  III 

where  pr  is  the  distance  from  the  image  of  the  line  source  to  the  receiver,  and  the 
+(-)  sign  is  used  for  the  case  of  a  magnetic  (electric)  line  source.  The  diffracted 
field  is  given  by  [4] 

m~jkp 

udif  (2.4) 

k  y/p 

where  uinc(Qe)  is  the  field  incident  on  the  edge,  and  the  diffraction  coefficient,  D», 
is  written  as 


D,  = - == 

A  2n\/2xk 


S[{ 


F\kVa-(4>  -  ^')]  J 

cotj^j^]  F(U"*a+(*+  *')]+ 

cot[id£±£l]  F(U-a-(^+^))  }J 
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The  transition  function  is  defined  by  [4] 


F(x)  =  2 jy/x  eJ*  /  e-jr2  dr 
J-Jz 


where  one  takes  the  principal  branch  of  the  square  root, 


(2.6) 


a±(P)  =  2  cos* 


2mrJV±  -  /? 


where  are  the  integers  which  most  nearly  satisfy  the  equations 


(2.7) 


2irniVi  -  0  —  ±ir 


(2.8) 


and 


n  =  2  -  WA/tt  (2.9) 

where  WA  is  the  wedge  angle.  Note  that  the  time  convention  is  assumed  and 
suppressed  throughout  this  study. 

The  subscripts  s  and  k  on  the  diffraction  coefficient,  D*,  correspond  to  the 

K 

cases  of  the  electric  and  magnetic  line  sources,  respectively  (these  arc  the  Boft 
and  hard  boundary  conditions).  The  four  terms  in  the  diffraction  coefficient  cor* 
respond  to  the  un  face”  incident  shadow  boundary  (ISB),  the  uo  face”  XSB,  the 
un  face”  reflection  shadow  boundary  (USB),  and  the  wo  face”  USB;  respectively. 
(The  o  and  n  faces  may  be  chosen  arbitrarily,  whereas  the  angles  <f>  and  4>'  must 
be  measured  from  the  o  face.)  The  transition  function,  as  introduced  by  Kouy- 
ournjian  and  Palhak  |4),  contains  a  Fresnel  integral  which  beastly  evaluated  using 
a  computer  algorithm.  The  magnitude  and  phase  of  the  transition  function  are 
shown  in  Figure  2.  Note  that  when  x  >  10,  F(x)  3  1.  The  region  near  a  shadow 
boundary,  where  F(x)  is  not  close  to  one,  is  called  a  transition  region.  In  a  tran¬ 
sition  region,  the  high-frequency  diffracted  field  can  not  be  considered  ray-optical. 


9 
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Figure  2:  Transition  function. 


The  L  parameters  are  known  as  the  distance  parameters,  and  are  given  by 


X4  = 
Lrn  - 
L”  = 


p'p 
p'  +  P 
P> 
Pe+P 
PcP 
Pc  +  P 


,  and 


(2.10) 

(2.11) 

(2.12) 


where  pj?  and  p”  are  the  caustic  distances  for  the  reflected  wavefronts  emanating 
from  the  edge  for  the  «  and  o  faces,  respectively.  (For  a  wedge  coni|M>sr  d  of  two  flat 
faces.  Pc  -  Pc  ~  p'-)  The  functions  o*(-)  arc  a  measure  of  the  angular  separation 
between  the  field  point  and  the  shadow  boundaries.  Finally,  note  that  the  following 
rules  must  be  observed: 


1.  Ll  >  A/10  (necessary  for  asymptotic  solution) 

2.  0<<M' <  n*  »  and 
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3.  If  (f)1  =  0  or  <j>'  =n7r,  one  has  grazing  incidence  and  must  multiply  the  diffrac¬ 
tion  coefficient  by  1/2. 

2.3  Three-Dimensional  Wedge  Diffraction 

In  the  previous  section,  the  scalar  fields  associated  with  a  cylindrical  wave 
incident  on  a  perfectly-conducting  two-dimensional  wedge  were  considered.  The 
problem  of  the  three-dimensional  wedge  (formed  by  two  plane  surfaces)  is  similar 
in  that  the  total  field  again  consists  of  incident,  reflected,  and  diffracted  fields. 
However,  the  three-dimensional  case  deals  with  a  point  source,  in  general,  and 
thus  illumination  by  an  obliquely  incident  spherical  wave.  Also,  the  incident  field 
may  be  arbitrarily  polarized.  To  simplify  the  expressions  for  the  reflected  and 
diffracted  fields,  two  ray-fixed  coordinate  systems  will  be  introduced. 

Consider  an  electric  field  E*  incident  on  a  perfectly-conducting  planar  surface, 

A 

as  shown  in  Figure  3.  Let  n  be  the  unit  normal  vector  to  the  surface,  /  be  the 
incident  unit  vector  from  the  source  to  a  reflection  point  Q%  on  the  surface,  and 

A  A 

R  be  the  reflection  unit  vector.  The  vectors  n  and  /  define  the  ordinary  plane  of 
incidence.  One  may  then  define  the  following  unit  vectors: 


„  n  x  J 

Til  =  - - 

|n  x  I\ 

(2.13) 

u[l  =  I  x  ux  ,  and 

(2.14) 

uj  =  R  x  ttx 

(2.15) 

where  1  and  ||  indicate  vectors  perpendicular  and  parallel,  respectively,  to  the 
ordinary  plane  of  incidence.  These  orthogonal  unit  vectors  form  a  ray*fixed  coor¬ 
dinate  system.  The  field  of  a  ray  traveling  in  the  incident  or  reflected  directions  is 
completely  specified  by  its  components  E±  and  £?||.  Thus,  the  reflected  field  E r(s) 
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at  a  distance  a  along  the  reflected  ray  from  the  reflection  point  Qr  may  be  written 
as 

£[j(*) 

E]_(a) 

where  E*(Qfl)  is  the  field  incident  at  Qr.  The  reflected  field  can  also  be  written 
as 

Er{s)*Ei{QR)>Re-jka  (2.17) 


1  0 
0  -1 


E^Qr) 

E[{Qr) 


(2.16) 


where  %  is  the  dyadic  reflection  coefficient  and  is  given  by 

U  =  ujufj  -  uitii  .  (2.18) 

To  examine  the  diffracted  field  in  this  problem,  consider  a  source  point  and 
observation  point  as  shown  in  Figure  4.  Let  the  distance  from  the  source  point  to 
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the  edge  diffraction  point  (Qe)  be  denoted  by  s',  and  the  distance  from  Qe  to  the 
observation  point  be  denoted  by  s.  Let  0'o  be  the  angle  of  incidence;  that  is,  the 
angle  formed  by  the  incident  ray  and  the  diffracting  edge.  When  an  incident  ray 
strikes  the  edge,  it  produces  a  cone  of  diffracted  rays.  The  axis  of  the  cone  is  the 
diffracting  edge,  while  the  cone  half-angle  (the  angle  of  diffraction),  /?0,  is  equal  to 
the  angle  of  incidence,  0'o.  The  equality  of  these  two  angles  uniquely  determines 
the  edge  diffraction  point,  Qg- 

Let  e  be  the  unit  vector  parallel  to  the  diffracting  edge,  I  be  the  incident  unit 
vector  from  the  source  to  the  diffraction  point,  Qe ,  and  D  be  the  diffraction  unit 
vector  from  Qe  to  the  observation  point.  The  plane  defined  by  e  and  I  is  the  edge- 
fixed  plane  of  incidence,  while  the  plane  defined  by  e  and  D  is  the  edge-fixed  plane 
of  diffraction,  as  shown  in  Figure  5.  One  may  define  the  following  unit  vectors: 


i<  -  -ix! 

|e  x  7| 

(2.19) 

0o  = 

(2.20) 

-a 

I 

X  X 

!l 

<-©. 

(2.21) 

0O  =  4>xb  . 

(2.22) 

The  unit  vectors  0'o  and  4>'  are  parallel  and  perpendicular,  respectively,  to  the  edge- 
fixed  plane  of  incidence,  and  form  the  ray-fixed  coordinate  system  fut  the  incident 
ray;  the  unit  vectors  Q0  and  ^  arc  parallel  and  perpendicular,  respectively,  to  the 
edge-fixed  plane  of  diffraction,  and  form  the  ray-fixed  coordinate  system  for  the 
diffracted  ray.  (Recall  that  in  the  case  of  2D  wedge  diffraction,  the  specification  of 
the  o  and  n  faces  was  arbitrary.  In  the  current  problem,  one  may  arbitrarily  assign 
the  direction  of  e.  Once  chosen,  however,  the  o  and  n  faces  are  also  determined. 
This  results  from  the  unit  vector  definitions  above,  and  the  conventions  shown  in 
Figure  4.) 
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Figure  5:  Ray-fixed  coordinate  system  (or  3D  diffraction. 


Using  these  vector  directions,  the  diffracted  field  at  the  observation  point  is 


written  as 

4s>) 

1 4(‘>  J 


~DS  0 

0  -Dh 


fytQs) 

tyiQs) 


P  e-jk. 


s(p  +  s) 


(2.23) 


where 


D,  = 


iw/4 


I-WSS-Al  M1^]  *“-*»-«]+ 

cot[ld^l]F(|.£0-W.f)]  J 

t{  cot  F(M,«+(*  +  f )]+ 

col[id^ill]F|Ha-w+W]  } 


(2.24) 


Since  the  wedge  is  formed  by  two  plane  surfaces,  the  distance  parameters  in  all 
four  terms  of  the  diffraction  coefficient  are  given  simply  by 


aa 


a  +  a 


7  sinz  (io  . 


(2.25) 


The  functions  a^(*)  and  F(>)  were  defined  earlier.  The  diffracted  field  may  also 
be  written  as 


Ed(s)~E  \Qe)'V, 


a(p  +  a) 

where  V  is  the  dyadic  diffraction  coefficient  and  is  given  by 


(2.26) 


V  =  -P’JoD,  -  i’iDk  . 


(2.27) 


The  geometrical  optics  and  edge  diffracted  fields  that  have  been  presented  are 
high-frequency  asymptotic  solutions  to  Maxwell’s  equations.  They  are  not  valid 
at  or  near  field  caustics  (although  one  may  allow  the  field  point  to  cross  a  caustic 
by  introducing  a  proper  phase  shift).  The  large  parameter  in  the  asymptotic 
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approximation  for  the  edge  diffracted  field  is  feL,  which  means  that  and  s' 

cannot  be  arbitrarily  small.  The  incident  field  is  assumed  to  be  slowly  varying 
at  the  edge  diffraction  point,  except  for  its  phase  variation  along  the  incident 
ray.  Also,  one  should  note  that  this  presentation  has  been  limited  to  a  wedge 
•  formed  by  two  plane  surfaces,  with  a  straight  edge.  Curvature  of  an  edge  primarily 

modifies  the  spread  factor,  while  surface  curvature  primarily  modifies  the  distance 
parameters.  A  much  more  detailed  and  general  presentation  of  this  subject  can  be 
found  in  [4]. 

2.4  Corner  Diffraction 

The  three-dimensional  edge  diffraction  solution  dealt  with  thus  far  has  been 
for  an  infinitely  long  edge.  In  an  actual  physical  problem,  the  scattering  object 
must  have  finite  length  edges.  Burnside,  Wang  and  Pelton  [5]  have  developed  a 
corner  diffraction  coefficient  which  compensates  for  the  termination  of  the  edge. 
Their  solution  is  based  on  the  equivalent  edge  currents  that  would  exist  in  the 
absence  of  the  corner.  They  employed  these  currents  in  the  radiation  integral, 
and  evaluated  the  integral  asymptotically.  This  analysis  was  characterized  by  a 
saddle  point  near  an  endpoint.  The  corner  diffraction  term  was  then  found  by 
appropriately  (but  empirically)  modifying  the  asymptotic  result.  Their  solution 
has  proven  to  accurately  predict  the  scattering  from  many  plate  geometries, 
s  The  corner  diffraction  geometry  is  illustrated  in  Figure  6,  where  a  corner 

exists  in  a  perfectly-conducting  fiat  plate.  The  corner  diffracted  fields  are  broken 
down  into  components  associated  with  each  individual  edge.  In  Figure  6,  the 
corner  diffracted  fields  associated  with  the  z  directed  edge  are  considered.  For  the 
source  and  observation  points  indicated,  there  is  an  edge  diffraction  point  indicated 
by  Q$.  As  the  observation  point  moves  in  the  -£  direction,  the  diffraction  point, 
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z 


Figure  6:  Corner  diffraction  geometry. 
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Qg,  also  moves  in  the  —  z  direction.  Eventually  it  moves  beyond  the  physical  limits 
of  the  edge,  and  the  edge  diffracted  field  stops  abruptly.  The  corner  diffracted 
field  compensates  for  this  discontinuity,  as  the  edge  diffracted  field  compensated 
for  the  discontinuities  in  the  GO  field.  Whether  Qg  is  on  the  edge  itself,  or  on  the 
imaginary  extension  of  this  edge,  the  corner  diffracted  field  will  exist. 

The  angles  /30,  /?c,  and  /?oc»  and  the  distances  a\  a",  ac,  and  a  are  illustrated 
in  Figure  6.  The  vector  directions  used  to  express  the  incident  and  diffracted  fields 
are  given  by 


—ex/ 

(2.28) 

j  = 

|e  x  i\ 

A  = 

$  x/ 

(2.29) 

ext) 

(2.30) 

i  = 

- - 7-  ,  and 

|e  x  D| 

Ac  = 

A  A 

fx9  . 

(2.31) 

The  corner  diffracted  fields  associated  with  one  corner  and  one  edge  of  a  perfectly- 
conducting  flat  plate  in  the  near  field  with  spherical  wave  incidence  are  given  by  [5] 


_  Em)  GaiQE)  I  e  y/sinftc  sin/jpc 

tfj,(Qc)  Ch(QE)  I  ( cos (3qc  -  cos/3c) 

F  (Mel* + A* -AM - 

& 


(2.32) 


where 


c.(Qe) 
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2\/2irk  sin  0O 


mAtzfi 


COJI 


F 

r 

kLca[*+Poc-Pc)\ 

F 

JcLca(*+(ioc-(k)\ 

(2.33) 
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Note  that  the  transition  function  was  defined  earlier,  while  the  distance  parameters 
are  expressed  as  follows: 


L  = 
Lc  = 


s' a" 


a*  +  s' 

scs 
sc  + a 


77  si«  0c 


(2.34) 

(2.35) 


and 


a(<j>  ±  <f>)  =  2  cos2  ^  . 


(2.36) 


The  corner  diffraction  coefficient,  C3(Qe)i  is  a  modified  version  of  the  edge  diffrac- 

h 

tion  coefficient,  D$(Qe).  The  modification  factor 
h 

La(<f>  ±  <ff)/ A 


[kLca(ir  +  0oc-  0c) 


is  a  heuristic  function  which  helps  to  avoid  the  abrupt  sign  changes  in  C3(Qe)  as 

h 

the  observer  passes  through  the  geometrical  optics  shadow  boundaries  of  the  edge. 

The  corner  diffracted  fields  associated  with  the  other  edge  are  treated  similarly. 
Thus,  depending  on  the  source  and  observer  locations,  there  may  be  edge  diffracted 
fields  from  either  or  both  edges,  but  there  will  always  be  a  corner  diffracted  field 
term  associated  with  each  edge.  The  total  effect  of  the  corner  i6  found  by  su¬ 
perimposing  the  contributions  associated  with  each  individual  edge.  This  topic  is 
discussed  further  and  applied  to  a  number  of  plate  geometries  by  Sikta  et  al  [10]. 

Consider  now  the  case  of  a  corner  formed  out  of  a  three-dimensional  wedge. 
The  corner  of  a  cube  is  an  example  of  this  type  of  geometry.  The  corner  diffracted 
fields  associated  with  one  corner  and  one  edge  of  such  a  perfectly-conducting  three- 
dimensional  wedge  in  the  near  field  with  spherical  wave  incidence  are  again  given 
by  Equation  (2.32),  where  the  comer  diffraction  coefficient  is  now  written  as 
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c . 

h 


x 


(2.37) 


— e~~J*/4 
2n\/2jrfc  sin  /30 
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'  -La* (4-4?) f\  ' 
kLca(x+(3oc-0c) 


+ 


cot 


La-(4>-4t)/X  1 

kLc<i(x  +0oc—@c) 


La+(*+4?)/ A  ' 

kLca(x+Poc—0c) 


+ 


cot 


kLca(x+l3oc-0c) 


The  corner  diffraction  solution  works  very  well  for  computing  backscattered 
fields.  In  [10],  the  authors  conclude  that  corner  diffraction  can  be  used  to  ob¬ 
tain  the  echo  area  of  rather  general  flat  plate  structures.  Unfortunately,  since  the 
corner  diffraction  solution  was  found  heuristically,  it  does  have  its  shortcomings. 
The  diffracted  fields  it  predicts  may  become  discontinuous,  particularly  for  some 
bistatic  geometries.  In  some  instances,  the  angles  0  and  <f>'  trigger  a  discontinuity 
that  the  modification  factors  cannot  entirely  eliminate.  These  values  of  <f>  and 
correspond  to  shadow  boundaries  of  the  GO  field.  In  terms  of  the  corner  diffrac¬ 
tion  solution,  these  are  false  shadow  boundaries,  since  the  corner  diffracted  field 
arises  to  smooth  the  edge  diffraction  shadow  boundary.  Nonetheless,  the  corner 
diffraction  solution  is  very  useful. 

Figures  8  to  10  demonstrate  both  the  validity  of  the  corner  diffraction  solution 
for  bistatic  scattering  cases  and  the  occurence  of  false  shadow  boundaries.  Consider 
a  plane  wave  incident  at  0*  =  45°,  <j>%  —  0°  on  a  2A  square  plate,  as  shown  in 
Figure  7.  Patterns  are  then  computed  for  fixed  4>a  and  for  0°  <  0s  <  360°.  Corner 
diffraction  results  are  shown  versus  moment  method  results,  both  are  calculated 
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in  dB  relative  to  a  square  wavelength.  The  actual  computations  were  done  by 
Brinkley  [11].  The  subscripts  on  cr  refer  to  the  source  polarization  and  receiver 
polarization,  in  that  order. 

In  Figure  8,  <f>a  =  0.1°,  which  places  the  observer  near  the  plane  of  incidence. 
The  nulls  in  the  moment  method  results  (for  the  99  polarization)  at  93  =  90°,  270° 
are  from  higher  order  effects  not  present  in  the  corner  diffraction  solution.  Agree¬ 
ment  between  the  two  solutions  is  excellent.  Note  that  there  is  a  slight  glitch  in 
the  corner  diffraction  results  at  6a  —  225°,  315°  which  is  due  to  difficulties  with  the 
modification  factor,  in  both  its  mathematical  expression  and  its  numerical  compu¬ 
tation.  This  difficulty  is  most  pronounced  in  the  plane  of  incidence  (<f>3  =  0°).  The 
value  of  0.1°  was  chosen  for  <f>9  to  lessen  the  dependence  of  the  displayed  results 
on  the  method  of  computation. 

Figures  9  and  10  depict  the  cross-polarized  and  co-polarized  results,  respec¬ 
tively,  for  <f>a  =  122°.  The  discontinuities  seen  in  the  corner  diffraction  solution 
are  due  to  false  shadow  boundaries.  Agreement  with  moment  method  results  has 
weakened  (in  comparison  with  the  case  where  the  observer  is  near  the  plane  of  inci¬ 
dence),  yet  the  dominant  scattering  behavior  is  predicted  by  the  corner  diffraction 
solution. 

The  diffraction  mechanisms  reviewed  here  serve  as  a  foundation  for  scattering 
from  dielectric  objects,  which  is  the  subject  of  the  next  chapter. 
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Figure  9:  Flat  plate  scattering  (cross-polarized  Helds),  <f> 3  -  122°,  (rom  (11). 
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CHAPTER  IH 


DIFFRACTION  FROM  DIELECTRIC  WEDGES  AND  CORNERS 

3.1  Introduction 

The  aim  of  this  chapter  is  to  develop  a  scattering  formulation,  based  on  the 
UTD,  for  a  dielectric  corner.  The  reflection  and  transmission  of  plane  waves  at  a 
dielectric-dielectric  interface  are  first  reviewed.  Next,  two  cases  in  which  dielectric 
modifications  have  previously  been  made  to  the  UTD  are  discussed.  These  were 
for  the  problems  of  edge  diffraction  from  a  two-dimensional  dielectric  slab  [6]  and 
from  a  two-dimensional  dielectric  wedge  [7].  By  extending  the  method  arrived  at 
in  [7|,  a  modified  UTD  solution  is  developed  for  three-dimensional  dielectric  edge 
and  corner  diffraction. 

3.2  Reflection  and  Transmission  at  a  Dielectric  Interface 

In  this  section,  attention  is  given  to  a  linearly  polarized  plane  wave  obliquely 
incident  on  a  planar  interface  between  two  media.  Medium  1  is  a  non-magnetic, 
lossless  dielectric;  medium  2  is  a  non-magnetic  and,  in  general,  lossy  dielectric. 

Consider  first  the  case  of  perpendicular  polarization,  where  the  electric  field 
direction  is  perpendicular  to  the  plane  of  incidence  (the  plane  formed  by  the  inci¬ 
dent  ray  and  the  interface  normal).  This  is  illustrated  in  Figure  11.  The  reflection 
coefficient,  JR±,  is  defined  as  the  ratio  of  the  reflected  electric  field  amplitude  to  the 
incident  electric  field  amplitude.  The  transmission  coefficient,  7^,  is  defined  as  the 


Figure  11:  Reflection  and  transmission  at  a  dielectric  interface,  perpendicular 

polarization.  . 


ratio  of  the  transmitted  electric  field  amplitude  to  the  incident  electric  field  ampli¬ 
tude.  These  Fresnel  coefficients  can  be  found  by  applying  the  boundary  conditions 
at  the  interface;  that  is,  by  enforcing  the  condition  that  the  tangential  electric  and 
magnetic  fields  be  continuous  at  the  interface.  Doing  so  results  in  the  following 

cos  -  sin2  9i 


R±  = 
T±  - 


cos 0*  -I-  ^(«2/ei)  -  sin20* 
2  cos 


and 


(3.1) 

(3.2) 


cos0*  +  7(c2/€i)  -  sin2  0* 

Note  that  Ti  =  1  +  R±. 

Consider  next  the  case  of  parallel  polarization,  as  shown  in  Figure  12.  The 
reflection  coefficient,  R^  is  defined  as  the  ratio  of  the  reflected  magnetic  field 
amplitude  to  the  incident  magnetic  field  amplitude.  The  transmission  coefficient, 
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> 


Figure  12:  Reflection  and  transmission  at  a  dielectric  interface,  parallel 

polarization. 


T||,  is  defined  as  the  ratio  of  the  transmitted  magnetic  field  amplitude  to  the 
incident  magnetic  field  amplitude.  Enforcing  boundary  conditions  as  before,  one 


finds  that 


(es/cpcoBfi*  -  yjjez/n)  "  8in20* 
(c2/«l)c°5 01  4  \/(«2/ei)  -  sin2 0‘ 

_  2(c2Al)cos  0* _ 

(62/«l)cos0*  4  \/(e 2A1)  ”  sin2  O' 


,  and 


(3.3) 

(3.4) 


Note  that  7y  =  1  4-  fly.  Since  the  incident  and  reflected  fields  are  in  the  same 


medium,  the  ratio  of  the  reflected  electric  field  amplitude  to  the  incident  electric 


field  amplitude  is  also  given  by  /?y. 

For  both  polarizations,  one  finds  that  0r  —  0*  and  sin0*  =  y/q/«2  sin0\  which 
are  known  as  Snell’s  laws  of  reflection  and  refraction,  respectively. 


29 


•  J»  .TJ.  Mi  Rj’  RM 


This  study  is  concerned  with  the  Fresnel  reflection  coefficients,  and  R^, 
for  reasons  that  will  become  apparent  later.  Note  that  R±  and  J2jj  are  complex 
for  complex  €3.  Also,  if  €2  is  real  but  less  than  ej,  and  R^  will  be  complex; 
however,  that  situation  will  not  arise  here.  Complex  €2  will  also  result  in  the  angle 
of  refraction  9t  being  complex.  (In  this  instance,  9t  is  clearly  not  the  true  angle 
of  refraction.  The  true  angle  of  refraction  can  be  found  in  a  manner  analogous 
to  that  done  by  Stratton  (12,  pages  500-505]  for  the  problem  of  refraction  in  a 
conducting  medium.  This  is  not  presently  of  interest.) 

Using  the  ray-fixed  coordinate  system  for  three-dimensional  reflection  shown 
in  Figure  3  (Section  2.3),  the  reflected  field  Er(a)  is  written  as 


£[j(») 

^Il 

0 

£f,(Qn) 

Ei(.) 

0 

_  Ei (Qr) 

(3.5) 


where  E*(Q#)  is  the  field  incident  at  a  reflection  point,  Q#,  and  a  is  the  distance 
along  the  reflected  ray  from  Qr. 

Figures  13  and  14  illustrate  the  behavior  of  |i2y  and  |JijJ  versus  angle  of 
incidence  for  three  values  of  C2»  Note  that  the  complex  permittivity  «2  can  be 
written  as  €2  =  -  jej},  where  €3  represents  the  lossy  nature  of  the  dielectric. 

When  the  material  is  lossless,  |ii|||  vanishes  at  the  Brewster  (or  polarizing)  angle. 
Also,  as  grazing  incidence  is  approached,  both  |j?|||  and  \R±\  approach  unity. 


3.3  Two-Dimensional  Dielectric  Slab 


The  first  application  of  the  UTD  to  a  dielectric  scatterer  was  developed  by 
Burnside  and  Burgener  [6],  who  studied  the  case  of  a  thin,  lossless  dielectric  slab. 
The  geometry  is  illustrated  in  Figure  15.  The  dielectric  slab  they  considered  was 
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Figure  13:  Plot  of  |J2jJ  vs.  6i  for  various  €2- 


Figure  14:  Plot  of  |Aj|]  vs.  0*  for  various  C2> 
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OBSERVATION 


Figure  15:  Geometry  for  diffraction  from  a  2D  semi-infinite  dielectric  slab. 
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of  any  thin  uniform  material,  but  subject  to  the  following  restrictions: 

•  The  diffractions  can  be  assumed  to  emanate  from  the  edge  point  Q%. 

•  Energy  leaving.  Qe  tangential  to  the  slab,  as  a  surface  wave  or  a  transmission 
through  the  endface,  can  be  assumed  negligible. 

•  Loss  within  the  slab  can  be  assumed  negligible. 


They  first  wrote  expressions  for  the  incident,  reflected,  and  transmitted  fields  of  a 
semi-infinite  dielectric  slab  (assuming  unit  magnitude  excitation)  as  follows: 

-•■~rrx-  in  reeionB  I  and  II,  and 

(3.6) 


u 


tnc  _ 


=  < 


ur*l  = 


0  in  region  III 

~7Z~ 


jjtotal  e  jn  reg|on  and 


0 


in  regions  II  and  III 


(3.7) 


and 


u 


iron  _ 


j 'total  e~*kpi 


(3.8) 


0  in  regions  I  and  II,  and 
--yfi-  in  region  III. 

Note  that  the  transmitted  (reflected)  field  is  written  as  the  product  of  a  total 
transmission  (reflection)  coefficient  and  the  field  that  would  have  been  incident 
upon  the  observer  from  the  source  (image)  in  the  absence  of  the  dielectric  slab. 
The  coefficients  Rioia*  and  Ttoi<^  account  for  the  multiple  interactions  that  take 
place  within  the  slab  to  produce  the  total  reflection  and  transmission  fields,  each  of 
which  is  an  infinite  sum  of  waves.  This  is  illustrated  in  Figure  16.  Rtotal  and 
are  each  expressed  as  an  infinite  series,  with  each  term  in  the  series  representing 
one  component  of  the  reflected  or  transmitted  field,  respectively.  /2*0*0*  and  Tioial 
are  found  by  proper  attention  to  the  various  phase  delays  and  coefficients  of  single 
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TOTAL  REFLECTED  FIELD 


TOTAL  TRANSMITTED  FIELD 


Figure  16:  Multiple  interactions  within  dielectric  slab. 


reflection  and  transmission  involved.  Note  that  plane  wave  propagation  was  as¬ 
sumed  for  the  field  incident  on  the  slab,  the  field  within  the  slab,  and  the  reflected 
and  transmitted  fields.  Thus,  neither  the  source  nor  the  observer  must  be  too  close 
to  the  slab. 

Now,  an  expression  for  the  diffracted  field  is  desired.  Comparing  the  dielec¬ 
tric  semi-infinite  slab  to  the  perfectly-conducting  case,  one  notes  that  the  discon¬ 
tinuities  at  the  two  shadow  boundaries  have  changed.  This  is  due  to  a  different 
reflected  field,  and  to  the  presence  of  a  transmitted  field.  Since  the  diffracted  field 
arises  to  smooth  out  the  discontinuities  at  the  GO  shadow  boundaries,  and  since 
these  discontinuities  have  been  scaled  by  and  1  -  T*0*0*  (  as  compared  to 
the  perfectly-conducting  case),  one  should  scale  the  diffracted  field  in  the  same 
manner.  This  was  the  approach  proposed  in  [6].  Recall  that  the  solution  for  the 
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(3.9) 


two-dimensional  perfectly-conducting  half-plane  is  given  by 

•  e~ikp 

ud‘t  ~u™(Qe)d/— 

where 

D$  =  D(<j>  -  <f>’)  T  D{4  +  4>')  (3.10) 

h  v  V  *  s— «v—~ ' 

ISB  term  RSB  term 

and 


. .  —  F\U4*±t) 

^  coe 


(3.11) 


The  diffraction  coefficient  is  then  scaled  to  [6] 


D,  =  (1  -  -  0#)  +  +  <f>')  .  (3.12) 

L  y  ii 

Note  that  the  solution  given  above  for  the  conducting  half-plane  is  a  special  case 
of  the  wedge  solution  given  in  Section  2.2,  found  by  setting  the  wedge  angle  to 

zero.  Also,  note  that  for  a  perfect  conductor  one  has  jTJ0*0*  =  0,  while  =  ^1. 

II  ii 

Thus,  the  proposed  dielectric  solution  reduces  to  the  correct  result  for  the  perfectly- 
conducting  case. 

The  validity  of  this  solution  was  verified  by  comparison  with  Moment  Method 
results,  for  the  case  of  Figure  17.  In  [6],  it  was  concluded  that  the  results  were 
accurate  for  p$  greater  than  a  wavelength  (which  stems  from  the  plane  wave  ap¬ 
proximation),  and  for  angles  of  incidence  <f>9  up  to  40°  (and  possibly  60°)  from 
normal  incidence  (to  avoid  exciting  a  surface  wave).  Finally,  the  theory  had  made 
the  assumption  that  the  thickness  D  was  less  than  1/10  of  a  wavelength.  How¬ 
ever,  based  on  comparison  with  measurements,  accurate  results  were  obtained  for 
thicknesses  up  to  1/2  wavelength. 


Figure  17:  Finite  2D  dielectric  slab  with  line  source  illumination. 

3.4  Two-Dimensional  Dielectric  Wedge 

In  the  previous  section,  a  dielectric  modification  was  made  to  the  UTD  en¬ 
abling  edge  diffraction  calculations  for  a  two-dimensional  dielectric  slab.  This 
section  considers  the  use  of  a  similar  modification  for  the  more  general  problem 
of  a  two-dimensional  dielectric  wedge.  This  problem  has  received  considerable  at¬ 
tention  by  a  number  of  authors.  However,  DeWitt[7]  reports  that  many  of  the 
published  solutions  to  this  problem  are  either  very  limited  or  too  complicated  to 
be  of  practical  use.  Alternatively,  RawlinB  [8]  has  developed  an  approximate  solu¬ 
tion  which  can  be  numerically  evaluated  fairly  easily.  He  derives  a  general  integral 
equation  for  an  infinite  dielectric  wedge  with  arbitrary  wedge  angle  and  permit¬ 
tivity.  A  Neumann  series  solution  of  the  integral  equation  is  then  found  through 
a  perturbation  technique.  Rawlins  takes  the  first  term  of  this  series  and,  using 
asymptotic  methods,  finds  explicit  expressions  for  the  diffracted  field  of  a  right 
angle  dielectric  wedge  under  these  constraints: 

•  The  incident  field  is  planar,  with  the  electric  field  vector  parallel  to  the  edge 
of  the  wedge. 


•  One  must  have  1  <  er  <  2  for  convergence  of  the  Neumann  series. 

o  The  incident  direction  is  in  the  quadrant  opposite  the  wedge. 

Rawlins’  solution  will  be  used  later  for  comparison  with  the  solution  to  be  presented 
here. 

The  dielectric  modification  to  the  UTD  discussed  now  was  proposed  by  De- 
Witt  [7],  as  an  extension  of  the  work  done  in  the  previous  section.  Since  it  is  in 
the  UTD  format,  it  will  be  numerically  efficient.  Furthermore,  as  a  simple  modi¬ 
fication  to  the  UTD,  it  will  not  have  the  constraints  mentioned  above.  The  first 
step  in  this  modification  is  to  consider  the  effect  that  the  dielectric  material  has 
on  the  discontinuities  at  the  shadow  boundaries.  It  is  immediately  apparent  that 
the  wedge  problem  is  much  more  complex  than  the  slab.  One  cannot  speak  of 
the  reflected  or  transmitted  fields  alone  in  that  the  multiple  interactions  which 
occur  within  the  wedge  are  of  a  more  complicated  nature.  However,  it  was  found 
in  (7]  that  the  following  method  works  very  well.  First,  due  to  the  vanishing  wedge 
thickness  at  the  edge,  energy  will  pass  through  the  tip  and  tend  to  smooth  the 
discontinuity  at  the  ISB;  thus,  the  1SB  terms  in  the  diffracted  field  are  considered 
negligible.  Second,  it  was  chosen  to  modify  the  RSB  terms  by  multiplying  them  by 
the  Fresnel  reflection  coefficient  for  the  initial  external  reflection  from  the  corre¬ 
sponding  wedge  face.  The  final  step  is  related  to  the  principle  of  reciprocity.  The 
solution,  as  presented  thus  far,  does  not  Batisfy  reciprocity.  This  is  true  since  the 
reflection  coefficient  depends  on  the  angle  of  incidence,  01,  which  changes  when 

source  and  observer  are  interchanged.  Thus,  in  [7],  the  Rx  were  calculated  by 

i 

replacing  O'  with  half  the  angle  between  the  incident  and  scatter  directions.  The 
resulting  formulation  for  the  field  diffracted  from  the  edge  of  a  dielectric  wedge  is 
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(3.13) 


h  y/p 

where  uinc(Qg)  is  the  field  incident  on  the  edge,  and  the  diffraction  coefficient, 
D» ,  is  written  as 


V01 

4-  Rx  cot 


F[fcIo+(^+^)] 


(3.14) 


^]F(Ma-(* +*')]}  • 


The  reflection  coefficients,  Rlt  are  given  by  Equations  (3.1)  and  (3.3),  respectively. 

H 

There  is  no  distinction  between  the  reflection  coefficient  for  the  o  and  n  face  terms, 
since  the  convention  used  for  the  angle  of  incidence  depends  only  on  the  incident 
and  scattered  directions. 

The  following  plots  compare  the  UTD  calculations  of  [7]  with  the  solution 
of  Rawlins  for  bistatic  scattering  from  a  two-dimensional  dielectric  wedge.  The 
electric  field  vector  of  the  incident  field  is  parallel  to  the  edge,  the  scattered  field 
is  computed  at  a  distance  of  50  wavelengths  from  the  edge,  and  the  relative  per¬ 
mittivity  of  the  wedge  is  equal  to  1.1.  This  value  of  er  is  chosen  because  Rawlins1 
solution  is  more  accurate  for  er  close  to  one.  The  three  figures  correspond  to  three 
incident  directions,  as  indicated.  Overall,  the  agreement  between  the  two  solutions 
is  very  good.  They  differ  most  significantly  along  the  wedge  faces.  Here,  one  would 
expect  Rawlins1  solution  to  be  more  accurate  since  it  uses  the  correct  boundary 
conditions.  On  the  other  hand,  the  modified  UTD  solution  may  be  less  accurate 
near  the  wedge  faces  since  the  presence  of  surface  waves  was  not  considered.  (Note 
that  later,  when  this  method  is  extended  to  pyramidal  absorber  scattering,  the 
incident  and  scatter  directions  never  lie  along  a  pyramidal  face.) 


Modified  UTD 


Modified  UTD 


Rawlins 


One  should  keep  in  mind  that  a  good  engineering  approximation  is  being 
sought  for  a  very  complex  problem.  The  method  proposed  in  [7]  has  proven  to 
work  very  well  for  the  cases  in  which  there  is  available  a  method  to  compare 
against.  This  modified  UTD  approach  will  serve  as  the  basis  for  the  work  in  the 
remainder  of  this  study. 

3.5  Three-Dimensional  Dielectric  Wedge 


In  this  section,  the  dielectric  modification  to  the  UTD  solution  proposed  in  [7] 
for  a  two-dimensional  wedge  will  be  extended  to  the  problem  of  an  infinite  dielec¬ 
tric  wedge  in  three-dimensional  space.  The  edge  diffracted  field  of  a  perfectly- 
conducting  three-dimensional  wedge,  as  depicted  in  Figure  4,  may  be  written  as 


Ed(a)  -  \AD{<j>  -  <j>')  +  BD{<j>  +  *')]  •  E HQb)J^ 


t — e^5 
+  *) 


where 


F[kLa+{<f>±  <*')] 


4-  cot 


F{kLa-(<t>±4>‘)] 


oi! 


and  the  dyadic  coefficients  A  and  B  are  given  by 


A  =  -hfa-ii' 

B  =  pj.-te'  . 


(3.15) 

(3.16) 

(3.17) 

(3.18) 


This  may  be  expressed  in  matrix  notation  as 


*tw 

Ep) 


E0'o(Qe) 

l4( Qb > 


*  s-** 


a(p  + j) 


(3.19) 
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where  A  and  B  are  matrices  given  by 
-1  0 

A  =  (3.20) 

0  -1 

1  0 

B  =  .  (3.21) 

0  -1 

The  goal  here  is  to  determine  the  dyadic  coefficients  A  and  B  that  will  properly 
scale  the  D(<f>  ^  <f>')  terms  when  the  wedge  is  dielectric.  The  approach  is  again 
similar  to  that  of  [6]  for  the  two-dimensional  dielectric  slab.  For  the  moment,  as¬ 
sume  that  R±  and  T±  represent  total  reflection  and  transmission  coefficients  for  the 

I!  II 

dielectric  wedge,  respectively.  The  discontinuities  at  the  reflection  shadow  bound¬ 
aries  will  now  be  investigated  to  determine  the  coefficient  B.  The  question  of  field 
polarizations  must  be  taken  into  consideration.  Recall  that  the  field  reflected  from 

a  planar  interface  is  most  conveniently  written  in  terms  of  its  ray-fixed  compo- 
i 

nents  u.l,  ujj,  as  was  done  in  Equation  (3.5).  On  the  other  hand,  the  diffracted 

A  .  ^  a 

field  formulation  deals  with  the  (30}  <f>  components  of  the  incident  field,  and  the 
$>»  4  components  of  the  diffracted  field.  To  find  the  discontinuities  in  the  reflected 
field  relative  to  these  field  polarizations,  consider  an  incident  ray  reflecting  off  a 
(planar)  wedge  face.  The  point  of  reflection  ib  infinitesimally  close  to  the  edge, 
so  that  the  reflected  ray  is  virtually  along  the  RSB,  on  the  lit  side  (as  shown  in 
Figure  21).  The  ordinary  plane  of  incidence  and  edge-fixed  planes  of  incidence 
and  reflection  are  indicated  in  Figure  22.  Define  a  to  be  the  angle  between  the 
ordinary  plane  of  incidence  and  the  edge-fixed  plane  of  incidence.  One  can  then 

A  .  v  . 

take  the  incident  field  expressed  in  its  ft  4>  components,  transform  them  into  u^, 

Uj|  components,  use  Equation  (3.5)  to  find  the  reflected  field  in  its  u^,  ujj  compo- 

*  • 

nents,  and  then  transform  this  reflected  field  into  its  0/>,  4>  components.  Doing  so, 
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Ordinary  Plane 
of  Incidence 


Edge-fixed  Plane 
of  Incidence 


Edge-fixed  Plane 
of  Reflection 


(View  along  reflected 
ray  toward  0R) 


Edga-fixad  Plane 
of  Reflection 


Edga-fixad  Plana 
of  Incidanca 


(View  along  Incident 
ray  toward  QM) 


Figure  22:  Pertinent  field  polarisation*  for  reflection  near  wedge  edge. 
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(3.22) 


one  finds  that  the  reflected  field  Er(s)  can  be  written  as 

E^s)  J2|j  cos2  a  —  i?j_  sin2  a  (fiy  +  R±)  sin  a  cos  a 

E^(s)  -( JZ||  +  ^l)  8^n  a  cos  a  — -^li  s*n2  a  +  COfl2  a 

x  [  1 
.  SjrW*)  . 

where  E %(Qr)  is  the  field  incident  at  the  reflection  point  Q#,  and  s  is  the  distance 
along  the  reflected  ray  from  Qr.  This  result  is  valid  for  any  incident  direction 
which  illuminates  the  wedge  face.  Since  the  reflected  field  is  zero  on  the  dark 
side  of  the  RSB,  the  2x2  matrix  in  Equation  (3.22)  represents  the  discontinuities 
associated  with  the  RSB  of  the  wedge  face  being  considered.  Thus,  the  scaling 
coefficient  for  the  D(<f>  +  <j>')  term  of  the  diffracted  field  is  written  in  matrix  form 


ii||  cos2  a  -  JR_l  8‘n2  a  (i?||  +  R± )  sin  a  cos  a 
B =  (3.23) 

-(i2j|  4- iij.)sinaco8a  sin2  ot  -+•  cob2  at  0 
where  the  subscript  indicates  that  the  values  vary  with  the  o  and  n  faces.  The 


angle  a  depends  on  the  o  and  n  faces  because  the  two  faces  have  different  normal 
vectors  and  thus  different  ordinary  planes  of  incidence.  Note  that  the  scaling  for 
the  o(n)  face  term  is  found  when  the  o(n)  face  is  illuminated,  but  will  be  used  for 
all  incident  directions  (regardless  of  whether  or  not  the  o(n)  face  is  illuminated). 
This  is  of  no  concern.  The  same  was  done  in  the  canonical  problem  of  the  two- 
dimensional  perfectly-conducting  wedge,  where  the  reflection  coefficients  are  ±1. 
When  a  certain  boundary  does  not  exist,  the  corresponding  term  in  the  diffraction 
formulation  is  typically  negligible. 

Following  a  similar  approach,  the  discontinuity  in  the  incident  shadow  bound¬ 
ary  (ISB)  terms  is  found  to  be  the  difference  between  the  field  on  the  lit  side  (the 
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incident  field)  and  the  field  on  the  shadow  side  (the  transmitted  field),  which  is 
expressed  as 


-1  +  T||  cos2  a  +  Ti  sin2  a  (Tjj  -  Ti)  sin  a  cos  a 
(Z]|  -  Ti)sinacosa  -1  4-T||8in2a  +  Tj.  cos2 a 


^{Qe) 

eUQe) 


e~jh$ 


(3.24) 


where  E*(Q£)  is  the  field  incident  at  Q$  and  s  is  the  distance  from  Q$  to  the 
field  point.  Thus,  the  scaling  coefficient  for  the  D((f>  —  <f>')  term  of  the  diffracted 
field  is  written  in  matrix  form  as 


A  = 


-1  +  T||  cos2  a  +  Tj_  Bin2  a  (Tj|  -  7j_)  sin  a  cos  a 
(7j|  —  IjJ  sin  a  cos  a  -1  +  T||  sin2  a  +  T±  cos2  a 


O 

n 


(3.25) 


where  the  subscript  again  indicates  that  the  values  vary  with  the  o  and  n  faces. 
This  procedure  was  first  used  in  [6)  for  the  problem  of  a  slab  in  three-dimensional 
6pace. 

Using  these  results,  and  rearranging  Equation  (3.19)  somewhat,  the  edge 
diffracted  field  of  a  three-dimensional  dielectric  wedge  can  be  written  as 


'  *&(•> 
.  E}(s) 


Da  Dh 
Dc  Dd 


&a(Qb) 

fyQs) 


s(p  +  a) 


where 


Da  «  (R\\  cos2  a  -  Rl  sin2  a)n  Dn(4>  +  4') 

+  (il||Cos2a  -  a)0D0((f>  +  <(f) 

4-  (-1  +  2’||  cos2  a  -f  sin2  a)n  Dn{<t>  -  ft) 
+  (-1  IT,,  cos2  a  +  T±  sin2  a)0  D0(<f>  -  <f>') 


(3.26) 


(3.27) 
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Db  = 


Dr.  = 


Dd  = 


[(12|j  +  R±)  sin  a  cos  a]n  Dn{<f>  +  <f>') 

+  [{R\\  +  R±)  sin  a  cos  a]0  D0(<f>  +  4>') 

+  l(r||  ~  Tl)  sin  a  cos  a]n  Dn{<l>  -  <f>') 

+  [(Tj|  -  TjJ  sin  a  cos  a]0  D0{<f>  —  <f>') 

-  [(i?||  +  -RjJ  sin  a  cos  a]„  Dn(<i>  +  <j>') 

-  [(iE||  -f  i?j_)sinacosa]0Po(^  + 

+  [(Tjj  -  T±)  sin  a  cos  o]»  !>*(#  -  <f>’) 

+  [(Tjj  -  Tjj  sin  a  cos  a]c  PD(^  -  <f>') 

(Rl  cos2  a  -  i2j|  sin2  a)n  Dn(^  +  <f> ') 

+  (il_LCos2a  -  i?j|  sin2  a)0  Po(<^  +  ^) 

+  (— 1  +  Tj|  sin2  a  -f  Tj_  cos2  a)n  Dn{<f>  —  <j>') 
+  (-1  +  Til  sin2  a  +  Ti  cos2a)0  D0{4>  ~  <t>') 


(3.28) 


(3.29) 


(3.30) 


and 


_e-j*/4 

DoW)  -  — 7== - 

w  2n\f2trk  sin  0O 


cot 


2  n 


F  (/?)]  . 


(3.31) 


Note  that  the  parallel  and  perpendicular  components  of  the  incident  and  diffracted 
field  are  cross-coupled.  This  is  due  to  the  dielectric  nature  of  the  wedge;  this  cross¬ 
coupling  does  not  exist  for  a  perfectly-conducting  wedge. 

The  preceding  formulation  is  included  for  the  sake  of  completeness;  it  aids  in 
visualizing  the  overall  effect  of  the  dielectric  (more  than  the  specialized  form  of 

this  solution  that  follows).  Recall  that  Rx  and  Tx  have  as  yet  not  been  defined. 

v  i 

Since  reflection  is  a  local  phenomenon  at  high  frequencies,  plane  wave  incidence 
is  assumed  and  Fresnel  coefficients  are  used.  The  final  step  here  will  be  to  use 
the  dielectric  modification  of  (7]  discussed  in  the  previous  section.  Thus,  Rx  will 

i 

consider  only  the  initial  external  reflections  off  the  wedge  faces.  Also,  recall  that 


the  usual  angle  of  incidence  is  replaced  by  half  the  angle  between  incident  and 
scattered  directions  in  the  computation  of  these  reflection  coefficients.  Therefore 

the  distinction  between  the  o  and  n  faces  for  #x  no  longer  exists.  (However,  the 

II 

distinction  between  the  o  and  n  faces  for  a  remains.)  The  Tx  are  set  to  unity,  which 

II 

is  equivalent  to  considering  the  ISB  terms  negligible.  With  these  conventions,  the 
edge  diffracted  fields  of  a  three-dimensional  dielectric  wedge  are  given  by 


Da  D}y 
Dc  Dd 


[  4«is)  j 


s(p+‘) 


where 


Da  =  (jR||  cos2  an  -  sin2  an)  Dn{<t>  4-  <t>') 

+  (#||  cos2  a0  -  R±  sin2  a0)  D0{<f>  4  <f>') 

Db  =  (#||  4#jL)sinancosanZ>n(^4^') 

-I-  (,R||  4  Ri )  sin  a©  cos  a0  D0(<j>  4-  <j>’) 

Dc  =  -  (#||  4-  #i)sinancosanDn(04-  <ff) 

-  (#||  4  R±)sina0cosa0Do(<l>  4-  ift) 

Dd  -  (R±  cos2  an  -  #||  sin2  an)  Dn{4>  4  tf) 

4  (#j_  cos2  a0  -  #||  sin2  a0)  D0(4>  4  ^) 

and  Do{4>  4  ^)  is  given  by  Equation  (3.31). 


(3.32) 


(3.33) 


(3.34) 


(3.35) 


(3.36) 


3.0  Dielectric  Corner  Diffraction 


At  this  point,  the  extension  to  dielectric  corner  diffraction  follows  almost 
intuitively.  The  corner  diffraction  solution  for  the  perfectly  conducting  three- 
dimensional  wedge,  presented  in  Section  2.4,  is  modified  here  just  as  the  edge 
diffraction  formulation  was  modified  in  the  previous  section.  Again,  the  dielectric 


conventions  of  [7]  are  used.  The  comer  diffracted  fields  associated  with  one  corner 
and  one  edge  of  a  finite  three-dimensional  dielectric  wedge  in  the  near  held  with 
spherical  wave  incidence  are  given  by 


’  ' 

rv 

-EJM 

Ca(QE)  Cb(QE ) 
Cc(Qe)  Cd(QE) 
Vsin/?csin/3oc 


fyLQc) 

[  fym 


o-J*/ 4 


(cos  0oc  -  cos  0C) 


F  [fcIco(ir  +  /3oc  -  0e)] 


y/2irk 

e-jks 


where 


and 


-c~J*/A 
2ny/2rrk  sin  0O 


(3.37) 


Ca(QE)  = 

(il||  cos2  On  -  R±  sin2  on)  Cn(<j>  4-  <f>') 

(3.38) 

4- 

(ilj|  cos2  a0  -  R±  sin2  a0)  C0{<f>  +  </>') 

Cb(QE)  = 

(i?||  4-  ill)  an  cos  an  Cn{4>  +  </>') 

(3.39) 

+ 

(Hjj  4-  ill)  8*n  °o  cos  00  C0(<j>  +  </)') 

Cc{Qe)  = 

— 

(iljj  4-  ili)smancosQtnCn(0  +  </>') 

(3.40) 

_ 

(il||  +  jRi)  sin  a©  cos  a0  C0{</>  4-  </>') 

Cd{QE )  * 

(ill cos2  “n  -  -Rjj  sin2  an)  Cn(</>  4-  4>') 

(3.41) 

4- 

(ill cos2  ao  ~  #||  Qo)  Co{<!>  4-  </>') 

(3.42) 


kLc<i{ir  4*  Poc  ”  0c)  j 
This  form  will  be  used  to  calculate  the  bistatic  scattering  from  the  tip  of 
an  absorber  pyramid.  It  is  applied  four  times,  once  for  each  of  the  four  edges 
associated  with  the  tip.  The  modeling  of  the  pyramid  and  specifics  on  how  this 
calculation  is  performed  are  the  subject  of  the  next  chapter. 
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CHAPTER  IV 


TIP  DIFFRACTION  FROM  PYRAMIDAL  ABSORBER 

4.1  Introduction 

In  this  chapter,  the  modified  UTD  solution  for  dielectric  corner  diffraction  will 
be  applied  to  an  absorber  pyramid.  The  modeling  of  the  pyramid  is  considered 
first,  and  the  calculation  of  several  variables  in  the  solution  is  discussed.  Results 
are  then  shown  for  various  monostatic  and  bistatic  geometries. 

4.2  Tip  Diffraction  Calculations 

The  pyramidal  tip  is  treated  as  a  corner  in  a  three-dimensional  wedge,  which 
has  four  planar  surfaces  and  four  edges  intersecting  at  a  common  point.  Recall 
that  the  corner  diffraction  solution  gives  the  fields  associated  with  one  corner  and 
one  edge;  thu6,  it  must  be  applied  to  the  pyramidal  tip  four  times. 

To  begin,  consider  the  conventions  illustrated  in  Figure  23.  The  tip  of  the 
pyramid  is  placed  at  the  origin  of  a  cartesian  coordinate  system.  This  will  serve 
as  the  primary  coordinate  system;  all  coordinate  systems  introduced  later  will  be 
defined  in  terms  of  this  primary  system.  The  pyramid  is  oriented  so  that  a  vector 
normal  to  face  1  or  3  would  have  no  y  component,  and  a  vector  normal  to  face  2 
or  4  would  have  no  £  component.  The  pyramid’s  shape  is  characterized  by  a 
single  angle,  a,  which  is  that  angle  between  any  two  adjacent  edges.  (Note  that 
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Edge  1  0 
Edge  4  N 

(FACE 


Face 

Face 

i) 


Figure  23;  Conventions  used  for  pyramid  geometry. 


(One  face  of  pyramid)  (Right  triangle  along 

pyramid  face) 

Figure  24:  Determination  of  the  pyramidal  wedge  angle,  WA. 

considering  the  tip  contribution  only  is  equivalent  to  the  pyramid  being  of  infinite 
extent.) 

First,  consider  the  wedge  angle,  WA,  a  variable  that  depends  only  on  the 
geometry  of  the  pyramid.  This  angle  is  most  easily  found  by  means  of  the  triangle 
shown  shaded  in  Figure  24.  This  triangle  has  two  vertices  at  opposite  corners  of 
the  pyramid  base,  while  its  third  vertex  lies  along  one  edge  of  the  pyramid.  This 
third  vertex  is  located  such  that  the  triangle  sides  adjacent  to  it  are  perpendicular 
to  the  pyramid  edge,  as  shown.  The  wedge  angle  is  the  angle  at  this  third  vertex  of 
the  triangle,  as  indicated.  Applying  the  law  of  cosines  to  one  face  of  the  pyramid 
yields 

b  =  yQ(l  -  cos  a)  (4.1) 

while  the  right  triangle  depicted  &lo»g  the  pyramid  face  yields 
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c  =  asina  . 


(4.2) 


Applying  next  the  law  of  sines  to  the  shaded  triangle  yields 
V2 b  c 

sin  WA  sin  7  (4.3) 

Substituting  Equations  (4.1),  (4.2),  and  sin  7  =  cos(WA/2)  into  Equation  (4.3) 
then  yields 

2a  yT  -cos  a  a  sin  a 

sin  WA  cos(WA/2)  *  (^*4) 

Since  sin  WA  =  2sin(WA/2)cos(WA/2),  one  finds  the  wedge  angle  to  be  given  by 

TIM  n  *  —  I  /  Vl  —  COSOi^ 

WA  =  2sm  1  ^  j  («) 

This  expression  is  of  course  valid  for  all  four  edges  of  the  pyramid. 

Next,  an  edge-fixed  coordinate  system  is  needed  for  each  pyramid  edge.  The 
coordinate  systems  chosen  are  rectangular;  each  has  an  edge  vector,  e,  directed 
along  the  edge  toward  the  pyramid  tip,  a  normal  vector,  n,  normal  to  the  o  face  of 
the  corresponding  edge,  and  a  bi-normal  vector,  e',  which  lies  along  the  o  face  and 
is  defined  by  e*  =  fi  x  e.  The  edge-fixed  coordinate  system  for  edge  1  is  depicted 
in  Figure  25.  Note  that  a  new  angle,  /?,  has  been  introduced,  which  is  the  angle 
formed  by  opposite  edges  of  the  pyramid;  it  will  be  useful  in  determining  the  edge- 
fixed  coordinate  systems.  Applying  the  law  of  cosines  to  each  of  the  two  triangles 
inset  in  Figure  25  and  equating  the  results  yields 


cos/3  «  2  cos  a  —  1  . 

Further  manipulation  reveals  that 

cos(/?/2)  =  ^/cos  a  ,  and 
sin(0/2)  v^2sin(a/2)  . 


(4.6) 


(4.7) 

(4.8) 


It  is  then  clear  that 


(4.9) 

(4.10) 


elz  =  cos(/?/2)  =  Vcosa  ,  and 
eh=ely  =  ~j*sin(0/2)  =  -mn(a/2)  . 

From  the  symmetry  present,  one  can  immediately  see  that  the  four  edge  vectors 


are  given  by 

ei  =  -*sin(a/2)  -ysin(a/2)  +  Zy/cosa  (4.11) 

e2  =  *8in(a/2)  -  ysin(a/2)  +  iv/cosa  (4-12) 

£3  =  x  sin(a/2)  +  y  sin(a/2)  +  zy/cos  a  ,  and  (4*13) 

I4  =  -x  sin(a/2)  +  y  sin(a/2)  4-  z \/cos  a  .  *  (4*14) 


Now,  let  us  consider  the  normal  vectors.  Since  nj  is  orthogonal  to  i\  and  to  the 
o  fexe  of  edge  1,  it  may  be  found  as 
yxei 


nl  = 


«  xA  +  zB 


(4.15) 


l»  *  ®ll 

where  this  computation  is  performed  numerically.  Note  that  A,B  >  0.  The 
remaining  normal  vectors,  due  to  the  symmetry  involved,  are  given  in  terms  of  A 
and  B  as 


ni  -  yA  +£2?  (4.16) 

n 3  -  '~&A  +  sB  ,  and  (4.17) 

«4  =  “ yA  +  sB  .  (4.18) 

Finally,  consider  the  bi-normal  vectors.  The  first  is  written  as 

c'j  =  n j  x  ej  -  xC  -  yD  -  zE  (4.19) 
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where  C,  D,  and  E  are  also  found  numerically.  Note  that  C,D,E  >  0.  Again, 
from  the  symmetry  present,  the  remaining  bi-normal  vectors  may  be  written  as 


e'2  =  xD  +  yC-zE 

ejj  =  -xC  +  yD  -  zE  ,  and 

£4  =  -xD-yC-zE  . 


(4.20) 

(4.21) 

(4.22) 


These  edge-fixed  coordinate  systems  will  be  used  later  to  determine  various  angles 
and  unit  vectors. 

Next,  consider  the  illumination  of  the  pyramid  by  a  spherical  wave  emanat¬ 
ing  from  a  source  point,  S.  Figure  26  illustrates  the  dielectric  corner  diffraction 
geometry  associated  with  the  tip  and  one  edge  of  the  pyramid.  The  source  and 
observation  points  are  in  the  far  zone  of  the  tip.  Thus  the  field  incident  at  the 
tip  is  approximately  planar.  Note  that  source  and  observer  are  assumed  to  recede 
into  the  far  zone  at  an  equal  rate;  that  is,  they  are  equidistant  from  the  tip.  Thus, 
while  the  location  of  Qe  >6  not  known,  there  is  no  ambiguity  in  the  value  of  0'o. 

The  computation  of  P'0  will  be  described  shortly.  Assume  that  the  incident  electric 

*  *  •  * 
field  is  known  in  terms  of  its  0,  <f>  components  and  its  incidence  angles  0\  <j>l\  as 

shown  in  Figure  27.  The  aim  is  to  find  the  scattered  field  in  a  direction  specified 

by  the  scattering  angles  0’,  4>a .  The  z,  y,  and  z  components  of  E*  are  found  by  a 

simple  transformation,  while  /  and  D  are  given  by 

/  s=  -zsin0‘cos^‘ -  ysiu0‘sin^  -  zcos0*  ,  and  (4.23) 

D  =  *  «n  0*  cos  +  y  sin  0*  sin  +  i  cos  0®  .  (4.24) 


Note  that  the  incident  and  scattered  directions  are  assumed  to  be  in  the  upper 
(2  >  0)  half-space. 
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To  implement  the  corner  diffraction  solution,  the  angles  <f>  and  4>'  must  be 
found  for  each  edge.  Consider  the  unit  vector  /  directed  toward  a  corner  diffraction 
point,  Qc ,  as  illustrated  in  Figure  28.  FVom  the  lengths  labeled,  it  is  clear  that  / 
is  expressed  in  the  edge-fixed  coordinate  system  as 


/  s  -e1  sin 0c  cos  -  nsin/3c  sin  4>  +  e  cos/3c  .  (4.25) 

In  a  similar  manner,  one  would  find  that 

D  -  e'sin/Joc cos 4>  +  nsin tocsin ^-fecos/?oc  •  (4.26) 

FVom  the  expression  given  for  /,  one  can  see  that  the  angle  of  the  m£/l  edge  is 
found  as 


#m  =  Uur' 


(427) 


i 
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(Unit  distance  from  Qc) 


/ 


Figure  28:  Determination  of  4>  and  $  at  an  arbitrary  edge. 


while  the  expression  given  for  D  is  used  to  find  <j>  of  the  edge  as 

= (1^)  •  <428> 
Note  that  the  minus  signs  in  the  numerator  and  denominator  of  Equation  (4.27) 
are  retained  to  allow  a  computer  algorithm  to  return  the  angle  in  the  proper 
quadrant.  By  definition,  0  <</>,<£'  <  2rr.  If  4>'m  or  <£*,  is  greater  than  wr,  then 
the  diffracted  fields  associated  with  the  edge  are  ignored,  since  that  edge 
is  shadowed  (the  incident  or  scattered  field  direction  is  inside  the  extension  of 
the  wedge).  Although  regions  of  space  do  exist  where  there  is  no  diffracted  field 
contribution  from  some  edge,  there  is  no  region  of  space  that  sees  no  contribution 
from  any  edge.  For  the  case  of  backscatter,  note  that  4hn  —  <t>m- 

The  reflection  coefficients,  Kj_  and  Ay,  are  found  by  Equations  (3.1)  and  (3.3), 
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respectively.  Recall  that  the  incidence  angle  used  in  these  calculations  is  half  the 
angle  between  incident  and  scattered  directions,  and  can  be  found  by 

0RC  =  ^cos~1(-/-I))  .  (4.29) 

These  reflection  coefficients  (R±  and  Ry)  apply  to  all  four  edges.  For  the  backscat- 
ter  case,  dfa  =  0. 

A.  A  A  A 

The  unit  vectors  <j> ,  /3Ci  4>,  and  /3oc,  which  were  defined  in  Section  2.4,  can 


now  be  more  conveniently  defined  (for  the  mifl  edge)  as 

<t>m  =  cos  <t>m  -  sin  4>m  (4.30) 

(4.31) 

<Am  =  nm  cos  <j>m  -  e'm  sin  <£m  »  and  (4.32) 

0ocm  —  *  &  •  (4.33) 

The  angles  0C  and  @oc  associated  with  the  mtli  edge  are  found  as 

Pcm  =  ,  and  (4.34) 

Pocm  —  cos  *(i)  •  em)  .  (4.35) 

For  backscatter,  0ocm  =  w  -  0cm- 


A  good  deal  more  effort  is  involved  in  finding  the  angle  0'o.  Consider  the 
geometry  illustrated  in  Figure  29.  R}  and  R*,  are  vectors  originating  at  the  tip 
and  terminating  at  the  source  and  observation  points,  respectively.  These  vectors 
must  be  of  equal  length,  although  their  actual  length  does  not  matter  here.  Thus 

*  A 

one  can  choose  Rj  ~  -J,  and  Ro  =  D.  Line  segments  are  drawn  from  the 
source  and  observation  points  perpendicular  to  the  edge  extension;  their  lengths 
are  labeled  ttm  and  1^,  and  the  points  at  winch  they  intersect  the  edge  extension 
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are  located  by  the  vectors  Rplm  and  Rp^,  respectively.  The  vectors  Rp4m  and 

Rd^_  can  be  found  as 
rom 


RP*m 

—  (R-s  *  fimjfim 

(4.36) 

RPcm 

=  (Ro  •  em)em 

(4.37) 

while  the  distances  km  and  lo™  can  be  found  as 


km  - 


Ra-RpJm|  »  and 
km  =  R-o  “  R-Pom 


(4.38) 

(4.39) 


Now,  consider  the  line  segments  drawn  from  the  source  and  observation  points 
to  the  edge  diffraction  point  QEm'  A  pair  of  similar  triangles  exists,  though  not 
coplanar,  as  indicated  in  the  inset  of  Figure  29.  Since 


3m  __  im  ~  3m 

km  km 

the  distance  sm  can  be  written  as 


»m  - 


kmkn 

km  +  km 


where 


(4.40) 


(4.41) 


Im  -  -  R,.,,,  |  =  (R.  -  R.)  •  im  .  (4.42) 

The  vector  R.^,  directed  from  the  tip  to  the  edge  diffraction  point,  is  then  clearly 
given  by 


R Dm  “  Rpfm  +  *mk 


and  the  unit  vector  is  given  by 


i0  s  Ji&L  “  R* 


(4.43) 


(4.44) 
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(4.45) 


With  lQEm  known,  for  the  mifl  edge  can  be  found  as 

P,om  =  co8~1(fQEm'im) 


and,  of  course,  (3om  =  Pom‘  In  the  case  of  backseat  ter,  this  reduces  simply  to 
0om=P’Om=  ^  ^  edges- 

Looking  back  to  the  dielectric  corner  diffraction  solution  presented  in  Sec¬ 
tion  3.6,  one  sees  that  the  angle  a  (the  angle  formed  by  the  ordinary  and  edge- 
fixed  planes  of  incidence,  not  to  be  confused  with  the  angle  o  of  the  pyramid)  is 
used  repeatedly  in  the  definition  of  C(i{Qe)i  Q>(Qu)>  Cc{Qe)i  and  ^(Qe)’  This 
angle,  although  useful  in  visualizing  the  diffraction  solution,  will  not  be  explicitly 
used  in  the  computation  of  the  corner  diffracted  fields.  Rather,  the  C{Qe)  may 
be  expressed  as 

CaiQs)  =  [tf M Rx  +  VW R)|] Cn{<t>  +  ^')  +  [ WaRi  +  WIR\\\ C0{<t>  +  V)  (4.46) 
Cb{QE)  =  [WJRi  -  ]<?«(*+*')  +  [VfWx"^§fil|Jc?o(0+0')  (4.47) 

Cc(Qe)  =  +  <t>‘)  +  [VWi  -  \W*ll]  W  -f  ^')  (4.48) 

and 

CdiQv)  *  [W5«x+ Cn(4  +  <t>')+  +  )£„(*+*')•  (4.49) 

The  variables  ^”,2^3,4  for  4he  vn^  edge  are  given  by 


V>lm  -  dlm  ■ 

V»L  “  ,and 

V’V  “  “im '  • 


(4.50) 

(4.51) 

(4.52) 

(4.53) 


The  unit  vector  u"  for  the  edge  is  given  by 


O 

Aft 

U±m  = 


Am  X  * 


A^x/ 


(4.54) 


Note  that  is  the  unit  vector  normal  to  the  o  or  n  face,  respectively,  of  the  mtli 
edge. 

Next,  consider  the  distance  parameters.  Recall  from  Section  2.4  that 


L  = 


Lr  = 


s's" 


y  4-  s' 


scs 

sc  +  s 


U  sin*2  /30  ,  and 


(4.55) 

(4.56) 


Thus,  L  and  Lc  increase  without  bound  as  source  and  observer  move  into  the 
far  zone.  For  this  reason,  the  transition  functions  F  (fclca( 7r  +  -  Pc)\  and 

F  |fcjLo±(0  +  <f>1) ]  can  be  6et  to  unity.  This  is  equivalent  to  saying  that  an  ob¬ 
server  in  the  far  zone  will  never  be  in  the  transition  regions  associated  with  the 
edge  diffraction  shadow  boundary  or  the  GO  reflection  shadow  boundaries,  respec¬ 
tively.  (This  is  due  to,  the  fact  that  the  transition  regions  collapse  to  the  shadow 
boundaries  in  the  far  zone;  an  observer  can  only  be  in  the  transition  region  by  be¬ 
ing  exactly  on  the  shadow  boundary.)  On  the  other  hand,  the  ratio  LjLc  remains 
finite  and  the  modification  factor  is  significant  in  the  solution.  Upon  considering 
Figures  26  and  29,  it  can  be  seen  that  LjLe  depends  not  on  the  absolute  values  of 
3,  3c i  s\  and  3M,  but  on  their  relative  values.  Thus  one  can  compute,  for  the 
edge, 


(L/Lc)i 


w  “  /  a  ' 


sm2/?om 


/  SCS  \ 

\ac  +  s) 


(4.57) 
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Note  that 


(4.58) 

(4.59) 

(4.60) 


sc  —  s  —  |Ra|  —  |Ro| 

Jm  =  |l<?£m|  =  |Ri>m-R-|  »  and 

3m  =  |D<3£J  =  * 

For  the  backscatter  case,  this  reduces  to  ( L/Lc)m  =  sin /?Cm. 

The  only  quantities  not  yet  mentioned  are  E'^Qc)  and  f^(Qc)’  These  are 
found  simply  by  taking  appropriate  dot  products. 

4.3  Sample  Calculations 

The  computations  discussed  in  the  previous  section  have  been  used  to  generate 
the  “Absorber  code”.  In  this  section,  the  Absorber  code  is  used  to  calculate  the 
scattering  from  a  single  pyramidal  tip  in  various  monostatic  and  bistatic  scattering 
cases.  The  conventions  defined  in  Figures  23  and  27  will  be  referred  to  here.  The 
following  figures  plot  the  ratio  of  the  amplitude  of  the  scattered  electric  field  to 
that  of  the  incident  electric  field.  Since  the  observer  is  assumed  to  be  in  the  far 
zone,  the  1  /a  term  in  the  corner  diffraction  formulation  is  ignored. 

First  consider  the  backscatter  case,  depicted  by  the  solid  curve  in  Figure  30, 
where  4>l  ~  -  0°  and  0°  <  0*  -  0*  <  90°.  The  incident  electric  field  is  0  po* 

lari  zed,  the  frequency  is  5  GHz,  and  the  pyrumidal  absorber  is  characterized  by 
o  a  21°  and  cr  =  1.45  -  j0.58.  As  0*  varies,  the  incident  ray  striking  the  tip 
sweeps  through  the  4>  -  0°  plane,  which  is  perpendicular  to  face  1 .  As  0*  increases 
through  10°,  a  small  discontinuity  is  seen;  here  face  3  drops  into  the  shadow  region, 
and  false  shadow  boundaries  (see  Section  2.4)  exist  for  the  o  face  of  edge  2  and 
the  n  face  of  edge  3.  As  0*  increases  through  79°,  the  source  is  broadside  to  face 
1,  and  a  singularity  is  Seen.  Note  that  a  singularity  in  these  far-zone  calculations 
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Figure  30:  Backscatter  in  the  <j>  =  0°  and  45°  planes  from  the  tip  of  an  absorber 

pyramid. 
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indicates  that  the  local  field  behavior  is  cylindrical  or  planar  (as  opposed  to  spher¬ 
ical).  (A  singularity  would  not  occur  in  a  near-zone  calculation,  where  all  distance 
variables  are  well-known,  all  transition  functions  are  included  in  the  calculation, 
and  the  1/s  spread  factor  is  not  dropped.)  In  the  solid  curve  of  Figure  30,  the 
singularity  occurs  when  the  source  (observer)  is  broadside  to  face  1.  At  broadside, 
the  observer  is  at  the  GO  RSBs  of  face  1  and  the  edge  diffraction  SBs  of  edges  1 
and  4.  The  incident  field,  assumed  planar,  produces  a  planar  reflected  field  from 
face  1.  To  smooth  the  discontinuity  caused  by  the  (planar)  reflected  field,  the  edge 
diffracted  field  must  behave  as  a  planar  field  near  the  RSB,  which  lies  along  the 
edge  diffraction  SB,  Thus,  the  corner  diffracted  field  must  behave  as  a  planar  field 
in  this  same  (backscatter)  direction  in  order  to  smooth  the  discontinuity  in  the 
edge  diffracted  field. 

The  second  (dashed)  curve  in  Figure  30  is  for  the  same  case,  except  that 
<f>x  =  =  45°.  At  Bx  —  0°,  both  curves  have  the  same  value  since  both  represent 

nose-on  incidence.  This  is  also  the  lowest  backscatter  value  on  the  plot,  indicating 
that  the  absorber  performs  best  at  nose-on  incidence,  as  expected.  In  the  second 
case,  as  6l  varies,  the  incident  ray  striking  the  tip  sweeps  through  the  <f>  =  45° 
plane,  which  contains  edge  1.  A  discontinuity  is  seen  at  6%  =  14°;  here  faces  3 
and  4  drop  into  the  shadow  region,  false  shadow  boundaries  exist  for  the  o  face  of 
edge  2  and  the  n  face  of  edge  4,  and  contributions  from  edge  3  are  ignored  since 
>  nir-  This  case  has  a  singularity  at  Bx  —  75°,  where  the  source  direction  is 
perpendicular  to  edge  1;  this  case  has  another  singularity  at  9X  =  90°,  where  the 
source  direction  is  perpendicular  to  both  edges  2  and  4.  In  these  instances,  the 
source  (observer)  is  along  edge  diffraction  SBs  only.  Edge  diffracted  rays  leaving 
the  edge  move  in  parallel  directions  to  the  (far  zone)  observer.  Thus,  the  edge  and 
corner  diffracted  fields  again  do  not  decay  as  a  spherical  wave,  and  a  singularity 
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results  in  the  calculations. 

Note  that  if  one  were  to  consider  a  finite  pyramid  (by  including  base  diffrac¬ 
tion),  these  singularities  would  not  exist.  Planar  GO  fields  would  emanate  from 
finite  faces,  and  planar  edge  diffracted  fields  would  emanate  from  finite  edges;  both 
would  radiate  at  one  precise  angle,  neither  are  seen  in  the  far  zone.  Of  course,  the 
tip  diffracted  term  does  not  change,  its  singular  component  cancels  with  a  like 
component  of  the  base  diffracted  field. 

Figure  31  illustrates  the  effect  of  the  pyramid  angle  a.  The  two  values  consid¬ 
ered  for  a  are  10°  and  21°;  other  parameters  are  indicated  in  the  figure.  Clearly 
the  tip  diffracted  fields  decrease  a  few  dB  with  decreasing  a.  Of  course,  for  con¬ 
stant  pyramid  height,  smaller  values  of  a  mean  more  pyramids  per  square  foot. 
Changing  a  also  affects  the  geometry  of  reflected  and  transmitted  fields. 

Next  the  effect  of  absorber  permittivity  is  illustrated  in  Figure  32,  where  the 
two  curves  correspond  to  er  =  1.45  -  j'0.58  and  er  =  1.38  -  j 0.2.  These  value* 
correspond  to  medium  and  light  doping,  respectively  [7].  Other  parameters  are 
indicated  in  the  figure.  Again,  the  tip  diffracted  fields  decrease  a  few  dB  with 
decreasing  permittivity  (note  that  it  is  the  imaginary  part  of  the  permittivity  that 
has  decreased  most  significantly).  Of  course,  attempting  to  use  this  lower  doped 
material  throughout  the  pyramid  to  reduce  the  scattered  fields  would  detract  from 
the  nature  of  the  absorbing  material.  However,  tapered  material  could  be  used  to 
solve  this  problem. 

The  frequency  under  consideration  has  no  effect  upon  the  nature  of  any  scat¬ 
tering  plot;  varying  frequency  can  only  shift  the  plot  up  or  down.  This  results 
from  the  pyramid  being  effectively  of  infinite  extent. 

A  bistatically  scattered  field  case  is  depicted  in  Figure  33,  where  the  inci¬ 
dent  direction  is  specified  by  $'  =  45°,  <f>1  =  0°  and  the  scattered  direction  by 
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Figure  31:  Backscatter  from  a  pyramidal  absorber  tip  in  the  <f>  -  0°  plane,  for 

two  values  of  a. 
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Figure  32:  Backscatter  from  a  pyramidal  absorber  tip  in  the  <f>  —  0°  plane,  for 

two  values  of  er. 
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0°  <  6s  <  90°,  <f>3  =  180°.  Again,  the  remaining  parameters  are  indicated  in  the 
figure.  The  first  discontinuity,  at  about  68  =  45°,  is  due  to  false  shadow  bound¬ 
aries  associated  with  the  o  face  of  edge  2  and  the  n  face  of  edge  3.  The  second 
discontinuity,  at  about  9s  —  55°,  is  due  to  false  shadow  boundaries  associated  with 
the  o  face  of  edge  4  and  the  n  face  of  edge  1. 

One  more  bistatic  case  is  illustrated  in  Figure  34,  where  0l  =  45°,  <f>1  =  0°, 
6s  =  45°,  and  0°  <  <f>a  <  360°.  Thus  the  incident  direction  is  fixed,  and  the 
scattered  direction  sweeps  through  a  conical  pattern  centered  on  the  z  axis.  There 
are  two  causes  for  the  various  discontinuities,  false  shadow  boundaries  and  edge 
shadowing  ( <f>  or  <j>'  not  in  the  range  from  0  to  nir  for  some  edge).  These  are 
specified  in  the  figure. 

While  the  discontinuities  seei  in  these  plots  may  seem  rather  severe,  one 
must  realize  that  the  levels  involved  are  extremely  low,  and  that  a  very  complex 
phenomenon  is  being  treated  approximately.  There  is  no  desire,  in  the  present 
study,  to  attempt  to  eliminate  these  discontinuities.  In  chap*er  VI,  comparisons 
made  between  calculations  and  measurements  for  both  a  single  pyramid  and  many 
pyramids  will  build  confidence  in  this  scattering  formulation. 
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Figure  33:  Bistatic  scattering  in  the  xr.  plane  from  the  tip  of  an  absorber 

pyramid 
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Figure  34:  Bistatic  scattering  along  the  9a  =  45°  cone  froir  the  tip  of  an  absorber 

pyramid. 


CHAPTER  V 


PYRAMIDAL  ABSORBER  SCATTERING  IN  A  COMPACT 

RANGE 


5.1  Introduction 

In  chapter  III,  the  formulation  for  dielectric  corner  diffraction  was  developed. 
Chapter  IV  then  described  how  to  apply  this  formulation  to  a  pyramidal  tip, 
thus  generating  the  single  pyramid  bistatic  scattering  code  (Absorber  code).  This 
chapter  investigates  the  scattering  from  the  pyramidal  absorber  material  in  the 
compact  range  chamber  that  occurs  during  normal  operation  of  the  range.  The 
method  of  calculation  is  discussed  first,  followed  by  sample  calculations. 

5.2  Method  of  Calculation 

The  scattering  into  the  target  zone  of  a  compact  range  from  the  pyramidal 
absorber  which  lines  the  ceiling,  side  walls,  and  floor  of  an  anechoic  chamber  is 
now  considered.  Figure  35  shows  typical  absorber  scattering  paths  for  ceiling  and 
floor  absorber.  Similar  paths  exist  for  absorber  along  the  side  walls.  The  analysis 
is  limited  to  tip  diffraction,  which  is  believed  to  be  the  dominant  scatterer  as 
the  incident  or  scatter  direction  moves  away  from  nose-on  incidence  [7,13].  This 
has  been  demonstrated  in  bistatic  measurements  where  scattering  was  measured  at 
various  angles  from  a  wall  of  pyramidal  absorber  [7].  In  any  case,  the  Absorber  code 
computes  only  tip  diffraction  and  will  be  used  here  to  predict  the  total  absorber 
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Figure  35:  Typical  absorber  scattering  paths  calculated  by  CHAMBER. 

scattering  that  takes  place  in  a  chamber. 

A  computer  program,  CHAMBER,  has  been  written  to  calculate  the  absorber 
scattering  in  a  compact  range  chamber.  CHAMBER  considers  a  simple  rectangular 
room,  as  depicted  in  Figure  36.  The  ceiling,  side  walls,  and  floor  are  lined  with 
pyramidal  absorber.  The  rear  wall  is  neglected  since  it  would  be  outside  the  range 
gate  in  a  well-designed  room.  A  focus-fed  semi-circular  parabolic  reflector  with  a 
rolled  edge  and  a  skirt  illuminates  the  room.  This  is  also  depicted  in  Figure  36. 
Note  that  the  origin  of  a  cartesian  coordinate  system  is  placed  at  the  reflector 
vertex.  This  is  henceforth  referred  to  as  the  “room  coordinate  system.”  A  few 
variables  used  in  CHAMBER  are  shown  in  Figure  36,  they  are  defined  as  follows: 

X-tips.ccil  =  Room  x*  coordinate  of  tips  of  ceiling  pyramids. 

X-tips  Jlor  =  Room  x1  coordinate  of  tips  of  floor  pyramids. 


76 


x_tipa_coil 


*_tips_?lor 


i_n(l,iu 


a  com 


b)  SIDE  VIEW 


Figure  36:  Auechoic  room  considered  in  absorber  scatter  calculations. 


y-tips 

z-refl.max 

z_room 


Room  y1  coordinate  of  tips  of  side  wall  pyramids. 
Maximum  z1  coordinate  of  range  reflector. 
Maximum  z1  coordinate  considered  in  the  room. 


It  is  through  these  variables  that  one  specifies  the  room  dimensions.  The  user 
must  also  specify  the  intended  location  of  the  target  zone  center,  along  with  ab¬ 
sorber  parameters,  feed  and  reflector  parameters,  and  timing  parameters.  This  is 
described  fully  in  Appendix  A,  which  explains  how  to  run  CHAMBER. 

CHAMBER  begins  by  choosing  patches  along  the  ceiling,  walls,  and  floor 
that  it  will  consider  for  its  scattering  calculations.  The  extent  of  tnese  patches 
roughly  coincides  with  the  variables  described  above.  However,  it  is  desirable  to 
have  patch  dimensions  that  are  multiples  of  0.5  ft.  Thus,  CHAMBER  finds  the 
following  variables: 


X„rOOt!)-!IUIl  = 

xj'oonunax  = 
y.rooin-min  = 
y_room~max  - 
z-rooin.iiiiii  - 
a -room-max  =3 


Smallest  multiple  of  -0.5  ft  >  x.tips  Jlor. 
Largest  multiple  of  0.5  ft  <  x-tips.ceil. 
Zero. 

Largest  multiple  of  0.5  ft  <  y-tips. 
Smallest  multiple  of  0.5  ft  >  z-reflanax. 
Largest  multiple  of  0.5  ft  <  ejroom. 


The  ceiling  patch  is  then  specified  by 

y_room.min  <  y1  <  y-rootn-inax 
z_room-inin  <  s'  <  z-room_max,  and 
z1  ~  x.tips.ceil 


(5.1) 

(5.2) 

(5.3) 
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while  the  side  wall  patch  is  specified  by 


x_room_min  <  x1 

<  x_room_max 

(5.4) 

z_room_min  <  z' 

<  z_room  .max,  and 

(5.5) 

* 

y' 

=  y-iips 

(5.6) 

<0 

and  the  floor  patch  is  specified  by 

A 

y_room_min  <  y1 

<  y_room_max 

(5.7) 

z-room-min  <  z1 

<  z_room.max,  and 

(5.8) 

J 

X 

=  x-tips_flor. 

(5.9) 

Note  that  these  patches  cover  half  of  the  ceiling,  one  side  wall,  and  halfol  the  fioor. 
However,  CHAMBER  does  effectively  consider  both  sides  of  the  room.  This  will 
become  clear  shortly. 

CHAMBER  then  divides  each  patch  into  “resolution”  cells  which  measure 
0.5*  x  0.5*.  Imaginary  grids  are  set  up  for  each  patch,  with  grid  points  at  the  center 
of  each  cell,  as  depicted  in  Figure  3?.  CHAMBER  assumes  that  a  pyramidal  tip  is 
present  at  each  grid  point,  and  all  calculations  will  be  made  with  this  imaginary  lip 
in  mind.  These  calculations  arc  then  assumed  valid  for  each  of  the  actual  pyramids 
within  the  corresponding  ceil. 

CHAMBER  now  determines  the  illumination  at  each  grid  point.  To  do  this, 
the  Semi-Circular  Compact  Range  Reflector  Code  (9)  has  been  taken  and  modified 
for  the  current  problem.  The  result,  subroutine  REEL-CODE,  provides  CHAM¬ 
BER  with  the  following  information  for  each  grid  point: 

•  The  (complex)  y\  and  z*  components  of  the  illuminating  GO  electric  field. 

(This  is  a  reflected  field  from  the  compact  range  reflector.) 
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Figure  37:  Imaginary  grid  along  ceiling  patch. 

•  The  corresponding  point  of  reflection  on  the  reflector. 

Note  that  the  Reflector  code  does  not  provide  accurate  phase  information  for 
the  fields  outside  the  target  zone,  due  to  an  approximate  method  used  in  finding 
reflection  points  along  the  rolled  edge  of  the  reflector.  (The  only  exception  to 
this  is  the  portion  of  the  floor  illuminated  by  the  skirt,  where  reflection  points  are 
easily  found.)  However,  phase  information  would  not  be  used  here  anyway,  due  to 
the  random  nature  of  the  tip  positions.  (Pyramidal  absorber  typically  comes  in 
2*  x  2*  blocks,  which  are  not  perfectly  aligned  when  placed  in  the  room.  Also,  the 
height  of  individual  pyramids  varies,  and  pyramids  tend  to  bend,  twist,  and  sag.) 
CHAMBER  does,  however,  use  the  phase  information  provided  to  determine  the 
relative  polarization  of  the  held  components,  which  is,  of  course,  essential.  Thus, 
at  this  point,  CHAMBER  has  the  following  information  at  cadi  grid  point: 
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•  The  amplitude  (and  sign)  of  the  x\  y',  and  z1  components  of  the  illuminating 
electric  field. 

•  The  incident  direction  of  the  field  illuminating  the  grid  point  (since  the  cor¬ 
responding  reflection  point  is  known). 

Next,  CHAMBER  must  compute  the  scattering  from  each  grid  point  to  the 
field  test  point.  The  choice  of  field  test  point  will  be  addressed  later  when  the 
various  options  available  in  CHAMBER  are  discussed.  For  the  moment,  assume 
that  the  test  point  is  fixed  somewhere  in  the  target  region.  The  Absorber  code  will 
be  used  to  compute  the  scattering  from  each  grid  point.  Recall  that  the  Absorber 
code  considers  a  pyramidal  tip  located  at  the  origin  of  a  cartesian  coordinate  sys¬ 
tem,  as  depicted  in  Figure  23.  This  is  referred  to  here  as  the  “pyramidal  coordinate 
system”.  The  Absorber  code  requires  the  input  variables  0\  £$,  E^,  63,  and  <j>3. 

These  pertain  to  the  pyramidal  coordinate  system.  Thus,  a  pyramidal  coordinate 
system  must  be  located  at  each  grid  point.  Note  that  the  z  axis  of  each  pyramidal 
coordinate  system  must  be  directed  into  the  room,  since  this  corresponds  to  the 
pyramid  being  directed  into  the  room.  Sample  celling,  wall,  and  floor  pyramid  co¬ 
ordinate  systems  are  shown  in  Figure  38,  along  with  the  room  coordinate  system. 

Consider  a  ceiling  pyramid  coordinate  system  whose  origin  is  at  some  obser¬ 
vation  point  (grid  point)  of  the  illuminating  field.  The  coordinates  of  this  origin 
in  the  room  coordinate  system  may  be  expressed  as 

(®V/,V)  a J*o,Sto,z<>)  .  (5.10) 

The  transformation  between  room  coordinates  (x',2/,2')  and  ceiling  pyramid  co- 
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Ceiling  Grid  Point 


Wall  Grid  Point 


A  K  y  "V 


■*  /■*./ 

-/y'yj 
/  x 


floor  Grid  Point 


Figure  38:  Typical  “pyramidal  coordinate  systems”  used  by  the  Absorber  Code. 


ordinates  (x,y>z)  is  given  by 


X  =  Zo  -  z 


y  -  yo-y1  *  w»d 


Z  =  Zo  -  X 


(5.11) 

(5.12) 

(5.13) 


Let  the  reflection  point  (on  the  reflector)  corresponding  to  this  grid  point  be  ex¬ 


pressed  in  the  room  coordinate  system  as 


(5.14) 


Then,  this  reflection  point  is  expressed  in  the  ceiling  pyramid  coordinate  system 


(x,y,z)  »  (zp  -  zr,y0  -  yr,x0-  xr)  . 


(5.15) 


Similarly,  if  a  grid  point  along  the  wall  is  specified  by 

(*\yV)  =  (*o,yo,*o)  (5.16) 

and  its  corresponding  reflection  point  is  specified  by 

(*\yV)  =  (*r>yr,2r)  (5.17) 

then  the  reflection  point  is  expressed  in  the  wall  pyramid  coordinate  system  as 

{x,y,z)-(z0~xr,Zo- Zr,yo~yr)  •  (5.18) 

The  analogous  result  for  a  reflection  point  which  corresponds  to  a  floor  grid  point 
is  expressed  in  the  floor  pyramid  coordinate  system  as 

(*,  y,  z)  -  ( yQ  -  yr,  z0  -  zr , -x0  +  *r)  .  (5.19) 

Thus,  every  grid  point  is  the  origin  of  a  pyramidal  coordinate  system,  and  the 
reflection  point  corresponding  to  each  grid  point  is  known  in  terms  >f  its  (®,y,e) 
coordinates  in  this  coordinate  system.  These  coordinates  can  be  used  to  find  the 
angles  of  incidence  required  by  the  Absorber  code  (0*, 

Next,  the  field  illuminating  each  grid  point  must  be  expressed  in  the  pyra¬ 
midal  coordinate  system  at  that  grid  point.  The  amplitude  (and  sign)  of  the 
illuminating  electric  field  is  known  in  room  coordinates  for  every  grid  point;  let 
these  (signed)  amplitudes  be  denoted  by  Exi ,  £y ,  and  Esi.  These  can  be  expressed 
in  the  pyramidal  coordinate  directions  at  a  ceiling  grid  poi<  t  through  the  following 
transformation 


Ex  a 

(5.20) 

Ey  » 

-Eyi  ,  and 

(5.21) 

E ,  = 

~EX>  . 

(5.22) 
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For  a  wall  grid  point,  the  transformation  is  given  by 

Ez  =  -Ex,  (5.23) 

Ey  =  -Ez<  ,  and  (5.24) 

Ez  =  -Eyi  (5.25) 

while  the  transformation  for  a  floor  grid  point  is  given  by 

Ex  =  -Ey,  (5.26) 

Ey  =  -Et,  ,  and  (5.27) 

Ez  =  Ex,  .  (5.28) 


Thus,  the  incident  field  components  (Eg,  E^)  required  by  the  Absorber  code  can 
be  found  for  all  grid  points  as 

E\  -  Ex  cos  0*  cos  ft  +  Ey  cos  O'  sin  ft  -  Ez  sin  O'  ,  and  (5.29) 

E\  =  -  Ex  sin  ft  +  Ey  cos  ft  (5.30) 

where  £*,  Eyy  Ez ,  0\  and  <j>'  pertain  to  that  grid  point. 

CHAMBER  finds  the  scattering  angles  (0a,  <f>a)  from  the  room  coordinates  of 
the  grid  point  and  field  test  point.  CHAMBER  then  provides  all  of  these  variables 
(01,  ft}  JEJ,  E^,  03 y  <f>3)  along  with  the  pyramid  angle,  a,  the  pyramid  relative 
permittivity,  «r,  and  the  frequency  under  consideration  to  the  Absorber  code  (for 
each  grid  point). 

The  Absorber  code  returns  the  complex  x,  y,  and  z  components  of  the  field 
scattered  from  every  grid  point  to  the  test  point  (with  4 he  1/s  spreading  factor 
being  ignored).  CHAMBER  uses  these  results,  along  with  the  distance  from  each 
grid  point  to  the  field  test  point  (to  account  for  spreading),  to  determine  the 
scattered  power  contribution  from  every  grid  point  to  the  field  test  point.  The 
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contribution  from  each  grid  point  is  weighted  by  the  number  of  pyramids  in  the 
corresponding  cell. 

This  explains  how  CHAMBER  computes  the  power  scattered  from  each  cell 
of  the  ceiling,  wall,  and  floor  patches  to  the  test  point.  Recall,  however,  that  these 
patches  lie  along  one  side  of  the  room,  and  it  was  stated  earlier  that  CHAMBER 
considers  both  sides  of  the  room.  To  understand  the  approach  used  to  account 
for  the  second  side  of  the  room  (y1  <  0),  consider  Figure  39.  Point  F  represents 
the  feed,  while  point  C  is  an  arbitrary  ceiling  grid  point,  and  R  is  the  reflection 
point  on  the  rolled  edge  corresponding  to  C.  Let  T  be  the  field  test  point  under 
consideration.  The  two  diagrams  in  Figure  39  correspond  to  two  locations  of  the 
test  point;  the  following  discussion  applies  to  both  cases. 

The  absorber  scattering  discussed  thus  far  corresponds  to  the  path  F-R-C-T. 
However,  for  every  ceiling  grid  point  C,  there  is  an  image  point  Cjm  from  which 
one  must  determine  the  power  scattered  to  the  test  point.  (Note  that  Cjm  has 
the  same  x\  s'  coordinates  as  C,  but  the  negative  of  the  y'  coordinate  of  C. 
RlM  and  Tjm  are  imaged  in  the  same  manner.)  This  second  scattering  path  of 
interest  is  F-Rim-Cjm-T.  Due  to  the  symmetry  present  in  this  problem,  the  power 
scattered  to  Tjm  via  path  F-R-C-Tjm  is  identical  to  the  power  scattered  to  T  via 
path  F-R|M-CiM-T.  Computing  the  scatter  via  the  latter  path  requires  less  of 
a  programming  effort  and  less  computer  run  time.  Thus,  CHAMBER  computes 
the  scattered  power  contribution  to  test  point  T  from  ceiling  point  C  via  path 
F-R-C-T,  and  computes  the  contribution  to  test  point  T  from  ceiling  point  Cjm 
via  path  F-R-C-Tim.  This  is  done  with  every  ceiling  grid  point,  60  that  the  entire 
image  of  the  ceiling  patch  16  considered.  The  wall  and  floor  patches  are  treated 
similarly.  (Note  that  when  the  test  point  moves  to  the  room  centerline,  the  test 
point  and  its  image  merge.  The  scattered  power  from  any  grid  point  and  its  image 
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igure  39:  Equivalent  absorber  scattering  paths  used  by  CHAMBER  for  the 

(y1  <  0)  side  of  the  room. 


* 


is  then  the  same,  and  need  be  computed  only  once.) 

Thus,  CHAMBER  computes  the  scattered  power  contribution  from  the  entire 
ceiling,  wall,  and  floor  patches  and  their  images.  However,  only  scattered  fields 
arriving  within  a  certain  time  window  are  considered.  This  is  handled  as  follows. 
It  is  assumed  that  the  chamber  operates  with  a  narrow  transmit  pulse  (in  the 
OSU  range  this  is  about  5  -  7  ns)  and  a  relatively  wide  receive  window  (in  the 
OSU  range  this  is  variable,  typically  set  to  20  ns).  The  target  zone  center,  the 
start  and  end  times  of  the  receive  window,  and  the  field  test  point  have  been 
specified  by  the  user.  The  time  reference  is  the  time  at  which  a  signal  would 
be  received  due  to  the  incident  plane  wave  scattering  from  a  point  target  at  the 
target  zone  center.  The  signal  scattered  from  some  grid  point  to  the  specified 
test  point  is  considered  only  if  that  signal  could  then  be  scattered  from  a  point 
target  at  the  test  point  and  arrive  at  the  feed  (via  the  reflector)  within  the  receive 
window.  This  is  illustrated  in  Figure  40.  The  reference  time  is  the  propagation  time 
along  path  F-R<j-TZC-R<j-F.  Assume  that  the  reference  time  is  70ns,  and  the 
specified  receive  window  start  and  end  times  are  -10ns  and  +10ns,  respectively. 
The  scattered  fields  from  the  grid  point  are  considered  only  if  the  propagation 
lime  along  path  F-Rjv-G-TP-Rp-F  is  between  60ns  and  80ns.  Other  timing 
approaches  arc  possible;  this  approach  computes  the  scattered  power  from  the 
absorber  which  could  be  scattered  from  some  feature  of  an  extended  target  (at 
the  specified  test  point)  and  arrive  within  the  receive  window,  thus  corrupting  the 
desired  return. 

Next,  there  is  the  question  of  addition  of  the  scattered  fields  from  the  various 
grid  points.  CHAMBER  computes  total  scattered  power  simply  by  adding  the 
scattered  power  contributions  from  the  individual  pyramids.  In  the  compact  range 
ceiling  scatter  experiment  described  in  the  next  chapter,  various  assumptions  of 
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Figure  40:  Timing  considerations  in  CHAMBER. 

field  addition  will  be  tested  in  the  calculated  curves. 

Finally,  note  that  all  scattered  power  calculations  are  normalized  to  the  power 
of  the  incident  plane  wave  in  the  target  zone.  Also,  if  one  requests  CHAMBER 
to  output  the  room  illumination  (from  the  range  reflector),  then  this  data  is  first 
normalized  in  the  same  manner. 

6.3  Sample  Calculations 

In  this  section,  CHAMBER  is  used  to  investigate  the  absorber  performance 
for  a  particular  compact  range  chamber  configuration.  As  required  by  CHAM¬ 
BER,  a  focus-fed  semi-circular  parabolic  reflector  with  a  rolled  edge  and  a  skirt 
is  considered.  The  focal  length  of  the  reflector  is  12  feet,  while  the  radius  of  its 
parabolic  section  is  11.5  feet.  The  roiled  edge  is  elliptical  with  a  major  axis  of 
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4  feet  and  a  minor  axis  of  1  foot.  The  feed  horn  is  tilted  above  the  horizontal 
by  25°.  The  tips  of  the  ceiling  pyramids  are  16.25  feet  above  the  vertex  of  the 
reflector  (at  x1  =  16.25  feet).  The  tips  of  the  side  wall  pyramids  are  20  feet  to 
each  side  of  the  vertex  of  the  reflector  (at  y1  =  ±20  feet).  The  tips  of  the  floor 
pyramids  are  3.5  feet  below  the  vertex  of  the  reflector  (at  x1  =  —3.5  feet).  The 
ceiling  and  side  wall  pyramids  are  characterized  by  a  =  24.5°,  er  =  1.45  -  j0.58, 
and  a  density  of  20.25  pyramids  per  square  foot.  The  floor  pyramids  are  character¬ 
ized  by  a  =  25.3°,  er  =  1.45  -  j0.58,  and  a  density  of  4  pyramids  per  square  foot. 
Note  that  the  height  of  the  skirt  is  not  specified,  but  it  is  assumed  to  extend  below 
the  tips  of  the  floor  pyramids.  The  target  zone  center  is  specified  to  be  on  the 
room  centerline,  6  feet  above  the  reflector  vertex, „  *  i  24  feet  downrange  from  the 
vertex.  In  the  room  coordinate  system,  this  is  ar)  =  (6,0,24).  The  relative 
start  and  end  times  of  the  receive  window  are  -10  ns  and  +10  ns,  respectively. 
The  chamber  is  assumed  to  operate  with  a  vertically-polarized  plane  wave  in  the 
target  region.  This  is  input  to  CHAMBER  by  specifying  the  appropriate  electric 
and  magnetic  dipoles  to  simulate  the  feed  (see  appendix  A). 

Now,  consider  a  lest  point  coincident  with  the  target  zone  center;  that  is, 
at  (zt,,y,,2>)  -  (6,0,24)  ft.  Figure  41  displays  the  total  scattered  power  versus 
frequency  at  this  test  point.  Also  shown  is  the  scattered  power  from  the  ceiling 
and  floor  (scattered  power  from  the  side  walls  does  not  arrive  during  the  specified 
time  window).  Clearly,  the  floor  is  the  dominant  scalterer  for  this  situation. 

Figures  42, 43,  and  44  display  the  variation  of  scattered  power  versus  test  point 
*\  y\  and  coordinate,  respectively.  In  all  three  figures,  the  frequency  is  fixed 
at  10  GHz.  In  Figure  42,  the  test  point  moves  upward  from  (x\yr,s*)  =  (2,0,24) 
to  (10,0,24)  ft.  The  floor  scatter  is  dominant  when  the  test  point  is  close  to  the 
floor,  while  the  ceiling  scatter  dominates  when  the  test  point  nears  the  ceiling. 
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ABSORBER  SCATTERED  POWER  (dB  rel.  to  incident 

plane  wave) 


Figure  41:  Absorber  scattered  power  vs.  frequency.  Test  point  at 
(*',y\a')  =  (6,0,24)  ft. 
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Owing  to  the  nature  of  the  Absorber  code,  some  floor  pyramids  yield  the  smalt 
jumps  present  in  the  range  2  <  x*  <  4  ft.  In  Figure  43,  the  test  point  moves 
across  the  room  from  (*',y',z')  =  (6, -5,24)  to  (6,5,24)  ft.  Clearly,  the  floor 
dominates  throughout  this  range.  Finally,  in  Figure  44,  the  test  point  moves 
downrange  from  (x',y',z')  =  (6,0,20)  to  (6,0,28)  ft.  Again,  the  floor  is  the  major 
contributor  to  the  absorber  scattered  power.  Note  that  increased  resolution  is 
provided  in  the  range  20  <  z1  <  21  ft.  The  singularity  at  20.3  feet  results  from 
Absorber  code  calculations  of  some  floor  pyramids.  Recall  that  the  Absorber  code 
is  designed  for  far-zone  calculations,  and  thus  sets  the  transition  functions  to  unity 
(see  Section  4.2). 

Now,  let  us  fix  the  frequency  at  10  GHz  and  the  test  point  at  the  target 
zone  center;  that  is,  at  (x',j /,*')  -  (6,0,24)  ft.  Figures  45,  46,  and  47  show  the 
absorber  illumination  as  a  function  of  position  along  the  ceiling,  side  wall  and  floor, 
respectively.  Only  half  of  the  ceiling,  one  side  wall,  and  half  of  the  floor  is  shown, 
since  the  illumination  is  symmetric  about  the  room  centerline.  To  understand 
the  discontinuities  present  in  these  figures,  one  must  realize  that  the  compact 
range  reflector  has  four  unique  sections.  Obviously,  there  is  the  parabolic  section 
which  produces  the  plane  wave  in  the  target  region.  The  skirt  is  a  cylindrical 
region  (independent  of  the  ar*  coordinate)  whose  y\  a'  behavior  matches  that  of 
the  parabolic  section  at  x'  -  0.  The  rolled  edge  consists  of  two  regions;  a  rolled 
edge  on  the  parabolic  section  (which  is  ^'-independent,  as  is  the  parabolic  section), 
and  a  rolled  edge  on  the  skirt  (which  is  x'-independent).  Figure  45  shows  that  the 
ceiling  illumination  is  most  intense  along  the  room  centerline  just  in  front  of  the 
reflector.  This  portion  of  the  ceiling  is  very  close  to  the  reflector,  so  that  not 
much  spreading  occurs  after  reflection.  There  are  no  discontinuities  in  this  plot 
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Figure  42:  Absorber  scattered  power  vs.  test  point  x'  coordinate.  Frequency 

10  GHz.  (y'  »  0,  z*  *  24  ft). 
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ABSOREER  SCATTERED  POWER  (dB  rel.  to  incident 
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Figure  44:  Absorber  scattered  power  vs.  test  point  z'  coordinate.  Frequency 

10  GHz.  ( x 1  =  6  ft,  y'  =  0). 
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since  all  reflection  takes  place  on  the  rolled  edge  of  the  parabolic  section.  On 
the  other  hand,  Figure  46  does  have  a  discontinuity  at  x1  =  0.  This  results  from 
the  reflection  point  having  moved  from  the  rolled  edge  of  the  skirt  to  the  rolled 
edge  of  the  parabolic  section,  thus  causing  a  sudden  change  in  the  spread  factor 
of  the  illuminating  field.  In  Figure  47,  the  discontinuity  at  y1  —  11.5  feet  is  due  to 
the  reflection  point  having  moved  from  the  skirt  to  the  skirt  rolled  edge.  Clearly, 
the  floor  is  strongly  illuminated  by  the  skirt.  The  absorber  scattering  to  the  test 
point  from  the  ceiling  and  floor  is  depicted  in  Figures  48  and  49,  respectively.  The 
absorber  scatter  is  symmetric  about  the  room  centerline,  because  the  test  point  is 
on  the  room  centerline.  Thus,  only  half  of  the  floor  and  ceiling  is  considered.  Now, 
one  can  see  that  the  receive  window  limits  the  portion  of  the  ceiling  and  floor  that 
contribute  to  the  total  absorber  scatter.  In  fact,  a  plot  of  absorber  scatter  from  the 
side  wall  is  not  included  since  none  arrives  during  the  specified  time  window.  The 
floor  scatter  shows  a  discontinuity  at  yf  =  11.5  feet,  due  to  the  discontinuity  in  the 
illumination  there.  These  absorber  scatter  plots  show  both  smooth  trends  in  the 
data  and  some  abrupt  changes.  This  is  due  largely  to  the  Absorber  code,  which 
predicts  such  behavior  as  pyramidal  incident  and  scatter  angles  vary.  Finally,  note 
that  CHAMBER  calculates  the  total  absorber  from  the  ceiling  and  floor  to  be 
-67.8  and  -55.8  dB,  respectively. 
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Figure  45:  Ceiling  illumination  V6.  ceiling  y'>  z'  coordinates.  (dB  rel.  to  incident 

plane  wave). 
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Figure  46:  Wall  illumination  vs.  wall  x\  z'  coordinates.  (dB  rel.  to  incident 

plane  wave). 
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Figure  47:  Floor  illumination  vs.  floor  y\  z1  coordinates.  (dB  rel.  to  incident 

plane  wave). 
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Figure  48:  Ceiling  scatter  vs.  ceiling  y',  z1  coordinates.  Test  point  at  (6,0,24). 

(dB  rel.  to  incident  plane  wave). 
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Figure  49:  Floor  scatter  vs.  floor  y',  2'  coordinates.  Test  point  at  (6,0,24). 
(dB  rel.  to  incident  plane  wave). 
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Now,  consider  a  configuration  identical  to  the  previous  one,  except  that  the 
y 1  coordinate  of  the  test  point  is  equal  to  5  feet.  The  test  point  is  thus  given  in 
room  coordinates  as  (x^y'^z1)  =  (6,5,24)  ft.  Figures  50  and  51  show  the  absorber 
scatter  to  this  test  point  from  the  two  sides  of  the  ceiling.  Figures  52  and  53  show 
the  absorber  scatter  to  this  test  point  from  the  two  sides  of  the  floor.  Again,  the 
side  walls  do  not  contribute  due  to  the  receive  window.  These  plots  do,  in  general, 
appear  as  one  would  expect  after  having  seen  the  results  for  the  previous  case. 
However,  the  lower  scatter  level  at  roughly  4  <  z*  <  12  ft,  3  <  y'  <  7  ft  was  not 
expected.  The  cause  was  found  to  be  due  to  the  pyramidal  scattering  angle  <f>3, 
which  slowly  passes  through  270°  in  this  region.  Edge  shadowing  and  RSB  effects 
cause  the  Absorber  code  to  predict  a  decreased  scattering  level  near  <f>3  =  270°. 
Finally,  note  that  CHAMBER  now  calculates  the  total  absorber  scatter  from  the 
ceiling  and  floor  to  be  -67.2  and  -56.0  dB,  respectively. 

Some  conclusions  may  be  drawn  at  this  point.  As  mentioned  earlier,  the  floor 
proved  to  be  the  dominant  scatterer  for  the  chamber  configuration  considered 
here.  Abo,  it  was  seen  that  the  walls  were  effectively  time-gated  out  in  all  the 
cases  considered.  In  Figure  44,  the  absorber  scatter  is  seen  to  be  fairly  constant 
in  the  range  20  <  z*  <  28  feet.  (Although  not  shown  here,  if  one  were  to  consider 
lower  values  of  z\  the  absorber  scatter  would  be  seen  to  increase  os  z1  decreases.) 
Most  importantly,  note  that  the  absorber  scatter  levels  predicted  here  are  much 
lower  than  typical  levels  of  diffraction  from  the  range  reflector  itself,  which  arc 
about  -30  dB  for  present  designs  [13]. 
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Figure  50:  Ceiling  scatter  vs.  ceiling  y\  z1  coordinates  (y1  >  0),  Test  point  at 
(6,5,24).  (dB  rel.  to  incident  plane  wave). 
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Figure  51:  Ceiling  scatter  vs.  ceiling  y',  z1  coordinates  (y*  <  0).  Test  poiut  at 
(6,5,24).  (dB  rel.  to  incident  plane  wave). 
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CHAPTER  VI 


EXPERIMENTAL  VERIFICATION  OF  PYRAMIDAL  ABSORBER 
SCATTERING  CALCULATIONS 

0.1  Introduction 

In  this  chapter,  absorber  scattering  calculations  are  verified  through  compar¬ 
ison  with  experimental  results.  Calculations  are  first  compared  with  backscatter 
measurements  of  a  single  pyramid.  Various  assumptions  on  the  addition  of  pyra¬ 
midal  fields  are  then  discussed;  these  will  apply  to  experiments  in  which  many 
pyramids  contribute  to  the  total  scatter.  Calculations  are  then  compared  against 
measurements  from  an  experiment  in  which  a  large  number  of  pyramids  were  scat¬ 
tering  under  the  same  bistatic  situation.  Finally,  calculations  are  compared  against 
experimental  data  in  which  the  bistatic  scattering  from  pyramidal  absorber  on  the 
ceiling  of  the  OSU  compact  range  was  measured. 

0.2  Backscatter  from  a  Single  Absorber  Pyramid 

Backscatter  measurements  of  a  single  pyramid  have  been  made  with  the  OSU 
compact  range  using  the  absorber  sample  depicted  in  Figure  54.  The  rear-facing 
pyramid  served  as  a  second  available  target,  and  as  a  means  of  balance.  Frequency 
scans  were  made  from  2  to  18  GHz,  in  10  MHz  increments,  for  three  orientations. 
Two  cases  are  presented  here;  these  are  4>  ~  0°,  8  —  0°  and  ^  =  0°,  9  a  45°.  The 
bandlimited  impulse  response  was  calculated  for  the  two  cases,  these  are  shown 
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Figure  54:  Absorber  sample  used  for  single  pyramid  backscatter  measurements. 
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Figure  55:  Bandlimited  impulse  response  of  single  pyramid;  4>  =  0°,  9  =  0°. 

in  Figures  55  and  56,  respectively.  The  expected  arrival  time  of  the  tip  return 
(-4.4  and  -3.1  ns  for  6  «  0°  and  45°,  respectively)  is  indicated,  along  with  other 
scattering  centers.  (Note  that  these  measurements  do  not  show  the  relative  levels 
of  the  tip  and  base  return  from  a  wall  of  pyramidal  absorber,  since  the  geometry 
of  a  pyramid  isolated  in  space  is  entirely  different  from  that  of  a  wall  of  pyramidal 
absorber.)  The  tip  return  is  larger  at  6  "  45°  than  at  $  -  0°,  which  was  predicted 
in  Figure  30. 

While  Figures  55  and  56  do  show  that  the  tip  response  begins  at  the  expected 
time,  it  is  diflicult  to  tell  when  it  ends.  For  this  reason,  the  impulse  response 
data  was  processed  further.  To  emphasize  signal  and  deemphasize  noise,  the  data 
was  squared  and  then  smoothed  with  a  moving  Hanning  window.  The  results 
are  presented  in  Figures  5?a  and  58a,  where  the  time  scale  is  now  from  -5  to 
-2  ns.  The  expected  time  of  the  tip  response  is  again  indicated  in  the  figures. 
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Figure  56:  Bandlimited  impulse  response  of  single  pyramid;  4>  -  0°,  9  -  45°. 

(Note  that  in  Figure  58a,  the  tip  response  appears  centered  at  -2.9  ns,  which 
corresponds  to  9  =  48°.  This  deviation  is  assumed  to  be  due  to  the  method  U6ed 
for  target  alignment,  which  was  believed  accurate  to  within  only  a  few  degrees. 
Of  course,  a  small  angular  error  would  noticeably  affect  the  tip  response  time  for 
9  ~  45°  ,  but  not  for  9  =  0°.)  Now,  one  can  visually  isolate  the  tip  return  and 
choose  appropriate  time  gate  parameters.  This  has  been  done  using  a  Kaiser-Bessc) 
window;  Figures  5?b  and  58b  show  the  gated  waveforms. 

These  plots  of  energy  temporal  density  were  considered  solely  to  find  ap¬ 
propriate  time  gate  parameters.  Going  back  to  our  calibrated  frequency  domain 
measurements,  one  can  perform  an  inverse  Fourier  transform  (where,  of  course, 
the  frequency  domain  data  is  not  windowed  before  the  1FT),  and  apply  the  proper 
time  gate  parameters  to  the  two  time  domain  signals.  Finally,  an  FFT  generates 
the  tip  response  in  the  frequency  domain  which  is  plotted  against  calculations  in 
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b)  Gated  Response 


Figure  58:  Energy  temporal  density  of  single  pyramid;  4>  =  0°,  6  =  45°. 
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Frequency  (GHz) 


Figure  59:  Calculated  vs.  measured  backscatter  from  a  single  pyramid;  <j>  —  0°, 

0  =  0°. 

Figures  59  and  60  for  9  =  0°  and  9  =  45°,  respectively.  In  each  figure,  vertical 
dashed  lines  indicate  the  range  of  valid  experimental  data,  due  to  the  processing 
performed.  The  measurements  show  that  the  scattered  signal  level  of  the  tip  is 
roughly  that  predicted  by  theory.  Better  agreement  can  not  be  expected,  since 
the  very  narrow  time  gates  used  limit  the  frequency  resolution;  in  addition,  the 
tip  return  is  a  very  weak  signal  being  measured  in  the  presence  of  much  stronger, 
signals. 

6.3  Comments  on  Addition  of  Pyramidal  Fields  and  Effective  RCS 

Up  to  this  point,  the  questions  of  how  the  fields  from  a  large  number  of  pyra¬ 
mids  add  to  produce  a  total  field,  or  of  how  an  “effective  RCS"  should  be  defined 
have  not  arisen.  In  Chapter  V,  it  was  decided  to  treat  the  pyramidal  absorber  in  the 


Frequency  (Gils) 

Figure  60:  Calculated  vs.  measured  backacatter  from  a  single  pyramid;  <j>  —  0°, 

$  =  45°. 

chamber  aa  an  incoherent  scatterer.  Thus,  the  values  of  power  found  from  various 
grid  points  were  simply  summed  (and  normalised)  to  produce  the  total  scattered 
power.  In  the  previous  section,  the  true  RCS  of  a  single  absorber  pyramid  was 
measured;  thus,  the  corresponding  calculations  were  straightforward.  However, 
the  experiments  described  in  the  following  two  sections  involve  the  scattering  from 
a  large  number  of  absorber  pyramids.  Neither  is  a  true  RCS  measurement.  Both 
involve  the  question  pf  how  the  fields  of  individual  pyramids  add  to  produce  the 
total  field.  For  these  reasons,  the  topics  of  pyramidal  field  addition  and  effective 
RCS  will  be  discussed  now. 

First,  lei  us  consider  a  few  cases  in  which  N  pyramids  are  cadi  illuminated  by 
the  same  field  E\  and  each  scatters  the  same  field  Btip  acot,  to  some  observation 
point.  Let  the  total  scattered  field  be  denoted  by  scat.  An  “effective  RCS” , 
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<reff,  is  then  defined  as 


CjjjtotaZ  scat  2\ 

4?rr  — E F —  J 


where  r  is  the  distance  from  the  tips  to  the  observation  point.  If  one  assumes  that 
the  N  fields  Etip  3Cat  add  in  phase  to  produce  Eiotal  3Cai,  that  is 

jgtotal  scat  _  jyjjjtip  scat  (6  2) 


then  it  follows  that 


<rcff  =  10  log  I  4tjt 


2  I  jgtip  scat 


+  201og  N  . 


On  the  other  hand,  if  one  assumes  that  the  N  fields  3cai  add  in  power  to 
produce  E<oio* acat,  that  is 

gtotoi  scat  2  _  jy  j  jj tip  »cot|2 
then  it  follows  that 


(ntip*cat  2\ 

4ffr2  — _ —  |  +  lOlogiV  . 


Note  that  the  assumption  of  field  addition  affects  only  the  log  N  term.  One  may 
want  to  base  the  assumption  of  field  addition  on  an  actual  measurement  of  pyra¬ 
mid  tip  position.  In  [7),  the  height  of  each  pyramid  on  an  absorber  block  of  36 
pyramids  was  measured.  Deviations  as  much  as  1  cm  were  observed.  Using  these 
measurements,  and  assuming  nose-on  plane  wave  incidence,  one  can  find  the  rel¬ 
ative  phased  the  tip  scattered  fields.  The  magnitude  of  the  sum  field,  $r»  was 
then  computed  as  follows  (7) 
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Figure  61:  Magnitude  (5r)  of  sum  of  fields  from  an  absorber  block  vs. 

frequency  [7]. 

where  <£j  is  the  height  of  the  itlv  pyramid,  k0  is  the  wavenumber  of  free  space,  and 
each  pyramidal  scattered  field  has  a  magnitude  of  unity.  The  results  (from  one 
block)  are  shown  in  Figure  61,  which  plots  the  magnitude  of  the  sum  field,  Sr, 
versus  frequency.  Phase  addition  corresponds  to  the  total  field  being  36  times  the 
field  of  a  single  pyramid  (5r  -  36),  while  power  addition  corresponds  to  the  total 
field  being  6  times  the  field  of  a  single  pyramid  (Sr  —  6).  At  the  lower  frequencies, 
the  curve  approaches  a  value  of  36,  as  tip  height  deviation  becomes  a  negligible 
fraction  of  a  wavelength.  Now,  this  data  will  be  used  to  appropriately  modify  the 
log  N  term  of  Equation  (6.3).  A  normalisation  constant  is  needed,  since  the  data 
in  Figure  61  was  found  for  a  block  of  36  pyramids.  The  resulting  expression  for 
effective  scattering  is  given  by 
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<reff  =  10  log  ^4?rr2 

(This  discussion  simply  serves  as  an  example.  If  one  were  to  take  this  approach, 
tip  position  measurements  of  the  material  in  question  should  be  made.) 

Now,  let  us  consider  several  blocks  of  absorber.  Again,  assume  that  each 
pyramid  is  illuminated  by  the  same  field  E*,  and  that  each  scatters  the  same  field, 
E tip  scat ^  some  observation  point.  The  total  scattered  field  produced  by  these 
pyramidal  fields  is  again  denoted  by  Etotai  3cat.  Since  absorber  blocks  tend  to  be 
misaligned  relative  to  each  other,  it  is  assumed  that  the  total  fields  of  the  absorber 
blocks  add  in  power  with  each  other.  However,  on  the  matter  of  the  addition  of 
pyramidal  fields  within  each  block,  three  assumptions  are  made: 

1.  Addition  in  power. 

2.  Addition  in  phase. 

3.  Addition  according  to  tip  height  measurements. 

Thus,  consider  M  blocks  of  absorber,  each  with  N  pyramids.  Under  the  first 

assumption,  one  finds  that 

C  fitipscol  2 
4*r 

while  under  the  second  assumption 

(I  ptip  scat  2\ 

4*r2 — —  I  -f  10log(A/JV2)  (6.9) 

and  under  the  third  assumption 

M/  =  10.og(4^|^n+10.og(MW{!^})  . 


j  -f  lOlog(MJV)  (6.8) 


E  tip  scat  2 

E* 


/log  Sr  \ 
\Iog  36  J 


log  AT 


(6.10) 


Note  that  Equations  (6.8)  and  (6.9)  are  special  cases  of  Equation  (6.10). 

Before  leaving  this  topic,  Equation  (6.10)  will  be  rewritten  for  the  more  general 
case  in  which  the  number  of  pyramids  may  vary  from  block  to  block.  Thus,  for 
M  blocks  of  absorber,  each  with  Ni  pyramids,  the  effective  RCS  (under  our  third 
assumption)  is  written  as 


<Teff  =  10  log 


jijttp  scat 


E* 


M 


+  10  log 


U=1 


(6.11) 


This  equation  applies  to  the  experiment  in  the  next  section. 

Note  that  a  major  restriction  to  all  of  the  cases  considered  in  this  section  was 
that  of  the  incident  and  scattered  fields  being  the  same  for  all  pyramids.  In  an 
experiment  which  measures  the  absorber  scattering  in  a  compact  range  chamber, 
this  will  clearly  not  be  the  case.  When  one  is  interested  in  computing  the  absorber 


scat 


term 


scatter  to  a  specific  test  point  (as  was  done  in  CHAMBER),  the  |Eioto* 
in  Equation  (6.1)  can  be  taken  to  represent  the  total  scattered  power.  However, 
since  each  pyramid  is  illuminated  differently,  there  is  no  obvious  interpretation  of 


the 


.m2 


E'l  term.  Fortunately,  one  does  not  encounter  this  difficulty  in  CHAMBER, 
since  only  the  scattered  power  computations  are  examined.  However,  the  calcula¬ 
tions  in  Section  6.5  force  one  to  consider  again  the  notion  of  an  effective  RCS.  This 
results  from  the  nature  of  the  experiment  considered  in  that  section,  however,  and 
is  best  explained  there. 


6.4  Bistatic  Scattering  from  a  Wall  of  Pyramidal  Absorber 


In  (7],  an  experiment  is  described  in  which  two  3'  diameter  parabolic  dish 
antennas  were  used  to  measure  the  bistatic  scattering  from  a  wall  of  pyramidal 
absorber.  The  wall  absorber  consisted  of  8**  tall  pyramids.  The  two  antennas  had 
broadband  TEM  horn  feeds  at  their  focii.  The  transmitter  produced  a  vertically- 
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polarized  electric  field.  The  experimental  setup  is  depicted  in  Figure  62.  The 
transmitting  antenna  was  positioned  at  an  angle  of  45°  relative  to  the  wall,  while 
the  receiver  was  positioned  at  a  variety  of  angles.  It  was  assumed  that  the  wall 
was  in  the  near  field  of  both  antennas.  Thus,  an  elliptical  region  (of  about  10 
v  sq  ft)  of  the  wall  is  assumed  to  be  illuminated  by  a  plane  wave.  The  receiver 

location  specified  by  <ps  =  45°  is  of  most  interest.  At  this  angle,  all  illuminated 
*  pyramids  have  (ideally)  identical  incidence  and  scattering  angles,  and  identical 

phase  paths  from  the  transmitting  horn  to  the  receiving  horn.  Thus,  depending  on 
how  much  variation  exists  in  tip  positions,  one  may  expect  the  scattered  fields  of 
the  pyramids  within  each  absorber  block  to  add  in  power,  in  phase,  or  somewhere 
in  between.  (Of  course,  the  pyramid  fields  could  interfere  with  one  another  and 
produce  a  total  field  less  than  that  predicted  by  an  adding  in  power  calculation. 
This  behavior  would  produce  nulls  at  certain  frequencies.  However,  the  power 
addition  is  considered  the  lower  limit  in  our  calculations.) 

Since  each  absorber  block  on  the  wall  measures  2'  x  2f,  a  different  number  of 
pyramids  are  illuminated  on  each  of  a  few  absorber  blocks.  The  three  assumptions 
of  field  addition  from  Section  6.3  will  be  used  for  pyramid  fields  within  each  ab¬ 
sorber  block,  while  the  total  scattered  fields  of  the  various  blocks  are  assumed  to 
add  in  power. 

Figure  63  shows  the  bandlimited  impulse  response  for  the  wall  scatter.  The 
.  location  of  the  expected  time  of  the  tip  and  base  returns  is  indicated.  Next, 

Figure  64  shows  the  corresponding  frequency  domain  data.  The  measured  data  is 
taken  from  [7],  while  the  calculations  were  made  using  the  Absorber  code  developed 
here.  Equation  (6.11)  is  used  to  find  the  three  calculated  curves;  this  is  done  by 
specifying  Sr  “  6  for  power  addition,  Sr  -  36  for  phase  addition,  and  using  Sr(f) 
from  Figure  61  for  addition  according  to  tip  height  deviation  (assumption  3).  The 
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WALL 


Figure  62:  Experimental  setup  of  bistatic  absorber  wall  scatter  measurement  (?]. 
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Figure  63:  Bandlimited  impulse  response  of  wall  scatter  measurement  (7]. 

third  assumption  agrees  best  with  the  measured  data.  Note  that  the  measurements 
include  both  tip  and  base  response.  Although  tip  diffraction  calculations  are  being 
plotted  against  measurements  of  both  tip  and  base  return,  recall  that  the  tip  return 
is  believed  to  be  the  dominant  scatterer.  Figure  63  supports  this  contention. 

0.6  Ceiling  Absorber  Scattering  in  a  Compact  Range  Chamber 

The  absorber  scattering  from  a  portion  of  the  ceiling  of  the  OSU  compact 
range  has  been  measured  by  means  of  a  large  flat  plate  placed  in  the  target  zone. 
This  work  was  the  effort  of  Young  and  Clerici  and  is  described  in  (14].  The  plate 
measures  6'  x  6'  and  can  be  rotated  in  azimuth  and  elevation.  The  experimental 
setup  is  depicted  in  Figure  65.  As  indicated  in  the  figure,  the  plate  is  assumed  to 
optically  reflect  the  incident  plane  wave  up  toward  the  ceiling.  (This  represents 
an  extreme  case  of  absorber  illumination  in  that  the  ceiling  absorber  just  above 


116 


&ef  f  (dB  rcl  to  eqnere  centimeter) 


Figure  64:  Comparison  of  frequency  domain  data  versus  calculations,  wall  scatter 
experiment  (experimental  data  from  (?)). 


Figure  65:  Experimental  setup  used  for  Plate  measurements. 


the  reflector  is  illuminated  by  the  incident  plane  wave.)  The  ceiling  pyramids 
then  scatter  in  all  directions.  Some  of  the  scattered  energy  will  strike  the  rolled 
edge  and  be  reflected  back  to  the  feed.  This  is  the  absorber  scatter  signal  under 
consideration. 

The  experiment  considered  here  used  a  plate  azimuth  angle  of  0°,  and  a  plate 
elevation  angle  of  73°.  This  corresponds  to  a  patch  of  ceiling  6  feet  wide  and  6.3 
feet  downrange.  (The  patch  is  actually  10.3  feet  in  extent  downrange,  but  the  last 
4  feet  are  lined  with  wedge  absorber  material.  The  wedge  material  is  not  believed 
to  scatter  significantly  in  the  direction  that  would  affect  tills  experiment.) 

Note  that  during  calibration,  the  range  normally  computes  the  following 


a  ss  lOlog 


(6.12) 


wh  .re  E*  represents  the  incident  plane  wave  and  E*  is  the  scattered  field  from 
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some  target  in  the  target  zone.  In  this  experiment,  the  range  will  receive  and 
process  our  ceiling  absorber  signal  as  though  it  had  originated  from  a  target  in  the 
target  zone.  Thus,  when  the  ceiling  absorber  signal  is  properly  time-gated,  the 
range  computes 

)  (6.13) 

/ 

where  E*  again  represents  the  incident  plane  wave,  and  |Ece^  ico*|2  represents  the 
total  scattered  power  of  the  ceiling  absorber  signal. 

Figures  66a  and  67a  show  the  bandlimited  impulse  response  for  the  cases  of 
vertical  and  horizontal  polarization,  respectively.  The  absorber  response,  based 
on  timing  calculations,  is  indicated.  A  Kaiser-Bessel  window  was  then  applied 
in  each  case,  Figures  66b  and  67b  show  the  gated  impulse  response  for  the  two 
cases.  Next,  an  FFT  was  performed.  Figure  68  shows  the  resulting  absorber 
response  in  the  frequency  domain,  along  with  three  calculated  curves.  In  these 
three  cases,  calculated  via  Equation  (6.13),  where  |EM^  "°*|2  is  found 

under  three  assumptions  of  held  addition.  The  lowest  of  the  calculated  curves  was 
found  assuming  the  fields  from  all  tips  added  in  power.  The  middle  curve  was 
calculated  assuming  the  fields  from  tips  within  one  square  foot  sections  added  in 
phase,  while  fields  among  sections  added  in  power.  The  upper  curve  was  calculated 
assuming  the  fields  from  tips  within  four  square  foot  sections  added  in  phase,  while 
fields  among  sections  added  in  power.  The  results  are  very  satisfying,  since  they 
roughly  predict  the  scattering  level  from  the  ceiling  patch,  and  since  they  follow 
the  frequency  behavior  of  the  measurements  as  well. 

(Note  that,  by  reciprocity,  equivalent  scattering  paths  exist  which  are  the 
reverse  of  those  shown  in  Figure  65.  These  were  the  paths  actually  considered 


Oceil  =  10  log  Uirr2 


E 


ceil  scat  I 


E* 
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when  the  ceiling  scatter  calculations  were  made.  This  was  done  by  a  computer 
code  similar  to  CHAMBER.) 
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IMPULSE  RESPONSE  IMPULSE  RESPONSE 
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a)  Total  Response 


b)  Gated  Response 

Figure  66:  Bandlimited  impulse  response  of  Plate  experiment,  vertical 

polarization. 
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b)  Gated  Response 

Figure  67:  Bandlimited  impulse  response  of  Plate  experiment,  horizontal 

polarization. 
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CHAPTER  VII 


SUMMARY  AND  CONCLUSIONS 

A  high  frequency  diffraction  solution  for  the  scattered  fields  from  a  pyramidal 
absorber  tip  has  been  developed  in  this  study.  This  solution  is  based  on  the 
perfectly-conducting  UTD  corner  diffraction  solution,  which  was  modified  to  treat 
the  case  of  a  corner  in  a  dielectric  three-dimensional  wedge.  The  modification 
is  modeled  after  that  done  in  [7]  for  the  two-dimensional  dielectric  wedge.  The 
resulting  formulation  was  then  used  to  generate  a  computer  code  for  the  calculation 
of  bistatic  scattering  from  an  absorber  pyramid.  The  validity  of  this  result  has  been 
verified  by  comparisons  with  backscatter  measurements  of  an  isolated  pyramidal 
absorber  tip.  Note  that  the  tip  term  was  extracted  from  the  measured  result 
by  transforming  the  frequency  domain  data  to  the  time  domain  and  gating  out 
the  appropriate  term.  The  calculated  and  measured  results  don’t  overlay,  but 
have  similar  levels.  One  should  not  expect  perfect  agreement  in  that  dielectric 
materials  have  many  potential  error  terms,  such  as  variation  of  dielectric  properties 
with  frequency  and  position.  Also,  the  tip  term  is  a  very  small  signal  which  was 
measured  in  the  presence  of  much  larger  signals.  Nevertheless,  the  measured  results 
do  tend  to  verify  the  theoretical  solution. 

The  validity  of  the  absorber  scattering  code  was  also  verified  by  comparisons 
with  bistatic  measurements  of  an  absorber  wall.  This  displayed  excellent  agreement 
between  calculation  and  measurements.  The  calculations  made  under  the  assump- 
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tions  of  power  and  phase  addition  (of  fields  from  individual  pyramids)  bounded 
the  measured  data  well.  The  calculated  curve  which  computed  pyramidal  field  ad¬ 
dition  based  on  tip  height  measurements  agreed  with  both  the  level  and  frequency 
dependence  of  the  measured  data. 

The  pyramidal  absorber  tip  solution  was  then  used  to  develop  an  anechoic 
chamber  analysis  for  compact  range  applications.  This  computer  code,  CHAM¬ 
BER,  determines  the  total  field  scattered  into  the  test  zone  by  the  absorber  lining 
the  chamber  walls,  ceiling,  and  floor.  (Note  that  CHAMBER  obtains  the  illumi¬ 
nation  of  the  absorber  from  a  reflector  analysis  code  which  treats  a  semi-circular 
reflector  with  a  rolled  edge  and  a  skirt.) 

A  specialized  version  of  CHAMBER  was  used  to  compute  absorber  scatter 
from  the  ceiling  of  the  Ohio  State  University  compact  range  chamber.  Compar¬ 
isons  of  these  calculations  against  measurements  from  the  plate  experiment  of 
Chapter  VI  illustrate  the  accuracy  of  our  predictions. 

CHAMBER  Has  been,  used  to  analyze  the  Ohio  State  University  compact 
range  chamber.  The  analysis  showed  that  the  dominant  scattering  came  from 
the  absorber  on  the  floor,  which  was  strongly  illuminated  by  the  skirt.  The  wall 
absorber  scattering  did  not  arrive  within  the  specified  receive  window,  for  cases 
where  the  test  point  was  along  ihr  room  centerline.  One  should  note  that  the 
fall  time  and  isolation  of  the  receive  window  deserve  careful  attention,  since  the 
scattered  power  (into  the  target  zone)  from  the  floor  and  ceiling  absorber  increases 
as  the  scatter  point  nears  the  back  of  the  room.  (This  can  be  seen  by  running 
CHAMBER  with  the  duration  of  the  receive  window  extended.)  Most  significantly, 
though,  the  absorber  scatter  levels  predicted  here  are  muck  tower  than  typical  levels 
of  diffraction  from  the  range  reflector  itself,  which  are  about  ~30  dB  (relative  to 
the  incident  plane  wave)  for  present  designs. 


With  this  code,  one  can  predict  the  pyramidal  absorber  performance  of  a 
compact  range  chamber  before  it  is  constructed.  In  fact,  the  various  hot  spots 
within  the  room  can  be  identified  so  that  appropriate  action  can  be  taken. 


APPENDIX  A 


PROGRAM  CHAMBER 


A.l  Introduction 

This  appendix  begins  with  a  description  of  hies  required  to  compile,  link,  and 
run  program  CHAMBER.  The  reader  is  then  stepped  through  the  operation  of 
CHAMBER,  and  the  output  generated  is  described.  Finally,  some  possible  error 
messages  will  be  discussed. 

A.2  Software  Requirements 

The  FORTRAN  code  for  program  CHAMBER  is  separated  into  two  hies, 
CHAMBER.FOR  and  CH-DECL.FOR.  The  former  contains  the  bulk  of  the  code; 
namely,  the  main  program  and  subroutines.  The  latter  simply  maintains  the 
variable  declaration  statements  and  common  blocks.  CHAMBER.FOR  inserts 
CH-DECL.FOR  wherever  necessary  with  the  “include"  command.  Thus,  the  hie 
CH-DECL.FOR  mu6t  be  present  when  compiling  CHAMBER.FOR. 

CHAMBER  also  uses  the  following  subroutines:  GEOMETRY,  BISTATIC, 
REFL-CODE,  and  PLOTTER.  Subroutines  GEOMETRY  and  BISTATIC  make 
up  what  has  previously  been  referred  to  as  the  Absorber  code.  Given  a  specified 
pyramid  angle,  a,  subroutine  GEOMETRY  computes  the  pyramid  wedge  angle, 
edge-fixed  coordinate  systems,  and  so  on.  It  passes  this  information  to  subroutine 
BISTATIC,  which  can  then  make  repeated  scattering  computations.  Of  course, 
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this  is  transparent  to  the  user  of  CHAMBER.  These  subroutines  are  both  in  file 
TIPBIST.FOR.  Next,  subroutine  REFL-CODE  is  the  modified  version  of  the  Semi- 
Circular  Compact  Range  Reflector  Code  [9]  that  has  been  mentioned  earlier.  The 
modifications  primarily  consist  of  conversion  from  a  program  to  a  subroutine,  and 
the  elimination  of  all  but  the  GO  field.  Also,  the  coordinates  of  every  observa¬ 
tion  point  and  corresponding  reflection  point  on  the  reflector  have  been  added 
to  the  output.  Again,  this  is  transparent  to  the  user.  File  REF3DT_SUB  con¬ 
tains  subroutine  REFL-CODE.  Subroutine  PLOTTER  creates  the  plots  output 

by  CHAMBER,  and  is  in  a  file  by  the  same  name. 

Finally,  when  running  CHAMBER,  two  data  files  must  be  present.  These 
are  CH_INPUT.DAT  and  FOR01T.DAT.  The  former  holds  the  input  data  that 
will  be  used  in  the  run  of  CHAMBER.  To  specify  the  input,  one  must  edit  this 
file.  This  is  more  convenient  than  entering  the  input  interactively,  especially  for 
multiple  runs  where  only  a  few  variables  change.  File  FOR017.DAT  contains  the 
cross-section  information  of  the  semi-circular  reflector,  which  is  used  by  subroutine 
REFL-CODE.  File  FOR017.DAT  is  created  by  program  SURFACE,  which  in 
turn  requires  file  FOR070.DAT  as  input.  (Note  that  some  data  specified  in  file 
FOR070.DAT  must  later  be  input  directly  to  the  reflector  code;  thus,  one  must  be 
careful  to  be  consistent.)  The  reader  is  referred  to  [9]  for  a  more  complete  discussion 
of  this  material;  however,  the  FOR070.DAT  file  that  created  the  FOR017.DAT  file 
used  in  the  sample  calculations  of  Chapter  V  is  as  follows: 
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Figure  69:  Range  reflector  cross-section  (y1  =  0),  skirt  not  shown. 
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0.,0. 

Figure  69  is  a  plot  of  the  data  in  FOR017.DAT;  note  that  this  represents  a  parabolic 
reflector  11. 5*  high  with  a  4'  x  l'  elliptic  rolled  edge  attached  at  the  top. 

In  summary,  the  following  steps  are  required  to  execute  program  CHAMBER: 

1.  Compilation  of  CHAMBER 

Files  required:  chamber.for,  ch.decl.for 
(VAX)  Command  :  fortran  chamber 

2.  Linking  of  CHAMBER 

Files  required:  chamber.obj,  tipbiat.obj,  refSdt  jtub.obj,  plotter.obj 
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(VAX)  Command  :  link  chamber,  tipbist,  ref3dt_8ub,  plotter 


3.  Running  of  CHAMBER 

Files  required:  chamber.exe,  chJnput.dat,  for017.dat 
(VAX)  Command  :  run  chamber 

With  the  executable  file,  CHAMBER.EXE,  one  can  begin  with  step  3. 

Note  that  some  options  of  CHAMBER  automatically  call  PLOTTER  to  plot 
various  output  data.  This  output  data  is  also  saved  in  output  files.  One  may  not 
want  CHAMBER  to  plot  the  results.  Also,  PLOTTER  is  likely  to  be  incompatible 
with  systems  other  than  that  which  it  was  written  on.  In  either  case,  the  calls 
to  subroutine  PLOTTER  can  easily  be  removed  from  CHAMBER,  and  no  results 
will  be  lost. 

A.3  Operation  of  CHAMBER 

In  this  section,  the  reader  is  brought  step-by-step  through  a  run  of  CHAM¬ 
BER.  The  CHJNPUT.DAT  file  that  was  used  in  this  run,  along  with  a  description 
of  each  input  variable,  is  as  follows  (note  that  all  linear  dimensions  are  in  feet): 

10.26  x'  coordinate  of  tipi  of  catling  pyramids. 

-3.5  x*  coordinate  of  tips  of  floor  pyranlds. 

20  y'  coordinate  of  tipa  of  aide  vail  pyramids. 

40  Maximum  z*  coordinate  of  intereat  in  the  room. 

8  x*  coordinate  of  target  zone  center. 

24  z'  coordinate  of  target  zone  center. 

24.5  Pyramidal  angle  (alpha)  -  ceiling  pyramids. 

20.25  lumber  of  pyramids  per  square  foot  (ceiling). 

1.46  heal  part  of  pyramid  relative  permittivity  (ceiling). 

-0.68  Xmag  part  of  pyramid  relative  permittivity  (ceiling). 

24.6  Pyramidal  angle  (alpha)  -  vail  pyramide. 

20.26  lumber  of  pyramide  per  square  foot  (vail). 

1.45  leal  part  of  pyramid  relative  permittivity  (vail). 

-0.68  Xmag  part  of  pyramid  relative  permittivity  (vail). 

26.3  Pyramidal  angle  (alpha)  -  floor  pyramids . 


lumbar  of  pyramids  per  square  foot  (floor). 

Real  part  of  pyramid  relative  permittivity  (floor) . 

Imag  part  of  pyramid  relative  permittivity  (floor). 

Focal  length  of  reflector. 

Radius  of  parabolic  section  of  reflector. 

Simulated  (1)  or  Measured  (0)  Feed 
Magnetic  dipole,  Xf  direction. 

Magnetic  dipole,  Yf  direction. 

Electric  dipole,  Xf  direction. 

Electric  dipole,  Yf  direction. 

Tilt  angle  of  feed  horn  (above  horizontal). 

Start  of  receive  vindov  (relative  to  time  reference). 

End  of  receive  vindov  (relative  to  time  reference). 

Most  of  these  are  self-explanatory,  and  have  also  been  mentioned  in  Chapter  V. 

Note  that  the  reflector  is  specified  as  having  a  simulated  feed.  The  following  four 

parameters,  described  as  magnetic  and  electric  dipoles,  are  used  to  simulate  the 

feed.  Details  on  this  are  provided  in  [9].  The  dipoles  chosen  here  (0, 1, 1,0)  specify 

vertical  polarization,  while  (1,0,0, -1)  specify  horizontal  polarization. 

The  run  begins  with  the  command  “run  chamber”.  CHAMBER  first  reads 
the  reflector  cros6-section  data  from  FOR017.DAT,  and  finds  the  maximum  x'  and 
z'  extent  of  the  reflector.  These  two  values  are  then  displayed,  along  with  resulting 
limitations  on  tip  position  coordinates,  as  follows: 

Physical  limit*  of  reflector: 
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Minimum  allovable  x  coord  of  ceiling  tip*....  =  15. Y  (ft) 
Minimum  allovable  y  coord  of  aide  vail  tip*..  °  15.7  (ft) 


Continue  (t/f)  7  T 
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Once  this  has  been  read,  and  appears  correct,  the  user  types  *t’  for  the  program  to 
continue.  CHAMBER  then  begins  reading  from  file  CHJNPCT.DAT.  The  room 
parameters  are  then  displayed,  along  with  information  on  the  resulting  ceiling, 
wall,  and  floor  patches,  as  follows: 

Room  parameters: 


X  coordinate  of  tipi  of  ceiling  pyramids....  =  16.25  (ft) 

X  coordinate  of  tips  of  floor  pyramids .  =  -3.50  (ft) 

Y  coordinate  of  tips  of  side  sail  pyramids..  =  20.00  (ft) 

Maximum  2  coordinate  of  interest .  =  40.00  (ft) 


The  patch  considered  along  the  ceiling  is  specified  by: 

Y  ranging  from  0.00  to  20.00  (ft) 

Z  ranging  from  4.00  to  40.00  (ft)  and  X  a  16.25  (ft). 

The  patch  considered  along  the  sail  is  specified  by: 

X  ranging  from  -3.50  to  16.00  (ft) 

Z  ranging  from  4.00  to  40.00  (ft)  and  Y  a  20.00  (ft). 

The  patch  considered  along  the  floor  is  specified  by: 

Y  ranging  from  0.00  to  20.00  (ft) 

Z  ranging  from  4.00  to  40.00  (ft)  and  X  a  -3.60  (ft). 


Continue  (t/f)  T  T 

Again  the  user  may  verify  the  data  and  type  V  for  CHAMBER  to  continue.  The 
target  zone  center  coordinates  are  then  displayed: 

Target  Zone  Center  Coordinates: 


x  a  6.00  (ft) 
y  =  0.00  (ft) 
x  a  24.00  (ft) 

Continue  (t/f)  T  T 
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Next,  the  absorber  pyramid  parameters  are  displayed: 

Pyramidal  absorber  parameters  (celling) : 


Pyramid  angle  (alpha) .  =  24.50  (degrees) 

lumber  of  pyramids  per  sq  It .  =20.25 


Relative  permittivity  ol  pyramid..  =  1.45  ~  j  0.5S 

Pyramidal  absorber  parameters  (aide  vail) : 


Pyramid  angle  (alpha) .  =  24.50  (degrees) 

lumber  oi  pyramids  per  aq  ft .  =20.25 

Relative  permittivity  of  pyramid..  =  1.45  -  j  0.68 


Pyramidal  absorber  parameters  (floor): 


Pyramid  angle  (alpha) .  =  25.30  (degrees) 

lumber  of  pyramids  per  sq  ft......  »  4.00 

Relative  permittivity  of  pyramid..  =  1.45  -  j  0.58 


Continue  (t/f)  ?  T 

Note  that  pyramidal  parameters  may  be  specified  separately  for  the  ceiling,  wall, 
and  floor  pyramids.  The  input  data  ends  with  the  reflcctor/feed  and  time  gate 
parameters,  as  shown  below: 

Reflector/Feed  parameters : 


Focal  length  of  reflector...... .  =  12.00  (ft) 

Radius  of  parebolie  section . .  »  11.60  (ft) 

Simulated  (1)  or  Measured  (0)  feed..  =  1 

Magnetic  dipole  (If) .  =  0.00 

Magnetic  dipole  (Yf) . . . =  1.00 

Electric  dipole  (If) . *  1.00 

Electric  dipole  (Yf) .  =  0.00 

Tilt  angle  of  feed . . .  «  26.00  (degrees) 


Time  gats  parameter* : 


Beginning  of  tine  gate....  =  -10.0  (na) 

End  of  tine  gate .  =  10.0  (na) 

Continue  (t/f)  ?  T 

Note  that  all  of  the  above  data  has  been  written  to  an  output  file,  CH_LOG.DAT, 

which  serves  as  a  log  of  the  entire  run. 

CHAMBER  now  computes  the  plane  wave  power  incident  at  the  target  zone 
center,  which  will  later  be  used  for  normalization  of  all  power  computations.  The 
round-trip  time  from  feed  to  target  zone  center  (via  reflector)  is  also  computed. 
A6  a  small  check,  CHAMBER  computes  the  plane  wave  power  incident  at  a  point 
2  feet  away  from  the  target  zone  center,  which  should  be  close  to  the  plane  wave 
power  at  the  target  zoue  center.  At  this  time,  CHAMBER  displays  the  following: 

Computing  poser  reference  (plese  save  poser)  end  tine 
reference  (round  trip  tiae  froa  teed  to  target  zone  center). 

Vote:  Poser  at  point  2  feet  eeey  froa  specified  target  zone  center 
ia  doen  by  0.11  dB. 

CHAMBER  is  now  ready  to  compute  the  illumination  of  the  ceiling,  walls,  and 
floor.  An  operating  frequency  has  not  been  specified  at  this  point,  but  is  irrelevant 
•inc«  we  do  not  consider  phase.  As  these  computations  are  made,  the  following 
messages  are  seen: 

Coaputlng  illuainatlon  of  calling... 

Coaputing  illuainatlon  of  aide  eall... 

Confuting  illuainatlon  of  floor... 


Coaputing  input  variable*  for  all  runa  of  subroutine  BISTATIC. . . 
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The  last  message  indicates  that  CHAMBER  is  converting  the  illumination  data 
for  each  grid  point  from  the  room  coordinate  system  to  the  pyramidal  coordinate 
system  at  that  grid  point. 

At  this  point,  CHAMBER  asks  the  user  to  specify  one  of  three  available  output 
options,  as  shown  below: 

AVAILABLE  OPTIONS: 


1)  Input: 
Output : 


2)  Input: 
Output : 

3)  Input: 
Output : 


Frequency . 

Test  point  in  target  region. 

Files  oi  ceiling,  vail  and  floor  illumination. 
FileB  of  ceiling,  vail  and  floor  scattering 
TO  the  test  point. 

Poser  (totals)  scattered  to  test  point. 

Frequency  range. 

Test  point  in  target  region. 

Poser  (totals)  scattered  to  test  point 
versus  frequency. 

Frequency . 

Range  of  test  points  in  target  region. 

Poser  (totals)  scattered  to  test  point 
versus  position. 


ENTER  the  desired  option  number:  1 

For  purposes  of  illustration,  option  1  has  been  chosen.  CHAMBER  then  requests 
the  user  to  specify  the  operating  frequency  and  test  point  location  (in  the  room 
coordinate  system).  These  questions  are  shown  (answered)  below: 

Enter  the  frequency  of  operation  (GHz) :  10 

Enter  the  Teet  point  X  coordinate  (ft):  6 

Enter  the  Test  point  Y  coordinate  (ft):  0 

Enter  the  Test  point  Z  coordinate  (ft):  24 
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Once  the  responses  have  been  made,  CHAMBER  lists  them  and  waits  for  verifi¬ 
cation  before  continuing,  as  shown  below; 

Frequency  of  operation  =  10.00  (GHz) 

Test  point  X  coordinate  ~  6.00  (ft) 

Test  point  Y  coordinate  =  0.00  (ft) 

Test  point  Z  coordinate  =  24.00  (ft) 

Input  OK  (t/f)  ?  T 

Note  that  the  following  is  added  to  CHJLOG.DAT: 

Option  1  was  chosen,  with  these  paraneters: 


Frequency  of  operation  a  10.00  (GHz) 

Test  point  X  coordinate  =  6.00  (ft) 

Test  point  Y  coordinate  =  0.00  (ft) 

Test  point  Z  coordinate  »  24.00  (ft) 

CHAMBER  next  computes  the  timing  data,  which  determines  whether  each 
individual  grid  point  is  included  in  the  scattered  power  calculations.  The  scattering 
from  the  appropriate  grid  points  is  then  computed.  As  this  is  done,  the  following 
messages  appear; 

Computing  timing  data. . . 

Computing  scatter  from  (entire)  ceiling  to  test  point... 

Computing  scattsr  from  (both)  side  sails  to  tsst  point... 

Computing  scatter  from  (entire)  floor  to  teet  point... 

Since  only  option  1  provides  the  room  illumination  as  output,  the  illumination 
data  is  now  normalised; 

formalizing  room  illumination  data... 
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The  results  of  the  scattered  power  computations  are  now  displayed  (and  written  to 
the  log  file).  The  contribution  from  the  side  walls  is  listed  as  —10,000  dB  because 
none  arrived  during  the  specified  receive  window.  CHAMBER  also  asks  the  user 
if  data  files  of  the  illumination  or  scattering  from  ceiling,  side  wall,  or  floor  are 
desired: 


Scattered  power  from  entire  ceiling..  =  -67.80  dB 

Scattered  power  from  side  walls .  =  -10000.00  dB 

Scattered  power  from  entire  floor....  =  -55.76  dB 

Total  scattered  poser .  =  -55.49  dB 


Enter  1  for  data  file  of  Ceiling  Illumination...  or 

Enter  2  for  data  file  of  Ceiling  Scatter .  or 

Enter  3  for  data  file  of  Hall  Illumination .  or 

Enter  4  for  data  file  of  Vail  Scatter..... .  or 

Enter  5  for  data  file  of  Floor  Illumination..,.,  or 

Enter  6  for  data  file  of  Floor  Scatter . .  or 

Enter  7  if  no  data  file  is  desired .  :  2 


Another  data  file  (t/f)  7  F 


To  illustrate  the  format  of  these  data  hies,  the  ceiling  scatter  file  has  been  written, 
and  the  beginning  of  this  file  is  shown  here: 

CEIL 

SCAT 


72 

40 

4.260000 

39.76000 

0.2600000 

19.76000 

-119.8876 

-119.8876 

-119.9283 

-119.9263 

-120.0026 

-120.0026 

-120.1166 

-120.1156 

-120.2626 

-120.2626 

-119.3672 

-119.3672 

-119.3183 

-119.3183 

-120.5702 

-120,5702 

-119.9973 

-119.9973 

-119.4919 

-119.4919 

-119.0696 

-119.0596 

-118.6570 

-119.9542 

-119.8186 

-119.6786 

-119.6563 

-119.4298 

-119.3396 

-119.2472 

-119.1503 

-119.0164 

-10000.00 

-10000.00 


-118.6579 

-119.0542 

-119.8186 

-119.8786 

-119.5563 

-119.4298 

-110.3396 

-119.2472 

-119.1603 

-119.0164 

-10000.00 

-10000.00 


The  labels  CEIL  and  SCAT  serve  as  identification  of  file  type.  The  numbers  72 
and  40  are  the  number  of  z1  values  and  y1  values  considered,  respectively.  The 
next  four  quantities  indicate  that  the  first  and  last  values  used  for  z'  are  4.25 
and  39.75,  respectively,  while  the  first  and  last  values  used  for  y1  are  0.25  and 
19.75,  respectively.  The  scattered  power  from  each  grid  point  is  then  written.  In 
column  1,  the  scattered  power  from  the  grid  point  whose  coordinates  are  specified 
by  the  first  z'  and  y1  values  is  first,  followed  by  the  grid  point  specified  by  the  first 
z 1  value  and  the  second  y1  value,  and  so  on.  When  the  last  y1  value  is  reached, 
the  second  z'  value  and  first  y1  value  is  considered.  The  second  column  holds  the 
scattered  power  contribution  from  the  image  of  the  corresponding  grid  point  in  the 
first  column.  In  the  case  shown,  the  two  columns  are  identical  since  the  test  point 
was  along  the  room  centerline.  Othei  data  files  have  the  same  format,  except  that 
illumination  files  will  have  only  one  column,  since  illumination  is  always  symmetric. 
Also,  for  wall  files,  the  coordinate  y'  is  replaced  by  x'.  These  data  files  can  be  used 
to  generate  any  type  of  plot  desired.  In  Chapter  V,  a  plotting  program  was  written 
to  generate  the  gray-scale  plots  shown  there. 
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CHAMBER  ends  the  run  with  this  reminder: 


YOU  SHOULD  PRIIT  THE  CH_LOG.DAT  FILE  !!! 


Now,  let  us  go  back  to  the  point  at  which  the  option  was  chosen,  and  select 
option  2.  Again,  CHAMBER  requests  user  input  at  this  point.  The  questions 
asked  are  shown  (answered)  below: 

EXTER  the  deaired  option  number:  2 

Enter  the  starting  frequency  (GHz):  2 
Enter  the  ending  frequency  (GHz):  18 
Enter  the  frequency  step  (GHz):  1 

Enter  the  Test  point  X  coordinate  (ft):  6 

Enter  the  Test  point  Y  coordinate  (ft):  0 

Enter  the  Teat  point  Z  coordinate  (ft):  24 


Once  the  responses  have  been  made,  CHAMBER  lists  them  and  waits  for  verifi¬ 
cation  before  continuing: 

The  start  frequency .  *  2.00 

The  end  frequency . =  18.00 

The  frequency  step .  a  1.00 

The  test  point  X  coordinate..  =  6.00 

the  teet  point  Y  coordinate..  «>  0,00 

The  test  point  Z  coordinate..  =  24.00 

Input  OK  (t/f)  ?  T 

Note  that  the  following  is  added  to  CH.LOG.DAT: 

Option  2  eas  chosen,  eith  these  parameters: 


(GHz) 

(GHz) 

(GHz) 

(ft) 

(ft) 

(ft) 


start  frequency . 

a 

2.00 

(GHz) 

end  frequency . . 

B 

18.00 

(GHz) 

frequency  step . 

e 

1.00 

(GHz) 

teat  point  X  coordinate. . 

5 

6.00 

(ft) 

teat  point  Y  coordinate. . 

a 

0.00 

(ft) 

teat  point  Z  coordinate. . 

a 

24.00 

(ft) 
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CHAMBER  will  first  determine  the  frequency  dependence  of  the  absorber 
scatter  of  the  ceiling,  side  walls,  floor,  and  of  the  room  as  a  whole.  This  is  done 
through  consideration  of  a  pyramid  with  parameters  identical  to  those  of  pyramids 
on  the  ceiling,  on  the  wall,  and  on  the  floor.  As  it  does  so,  this  message  is  seen: 

Computing  frequency  dependence  of  room  scatter . . . 

Next,  CHAMBER  computes  the  timing  data  (which  is  valid  for  all  frequencies) 
and  the  absorber  scatter  at  the  starting  frequency.  As  this  is  done,  the  following 
messages  are  displayed: 

Computing  timing  data... 


Computing  scatter  from  (entire)  ceiling  to  test  point... 
Computing  scatter  from  (both)  side  sails  to  test  point... 
Computing  scatter  from  (entire)  floor  to  test  point... 


CHAMBER  then  uses  the  frequency  dependence  information  to  determine  the 

absorber  scatter  at  each  frequency  specified.  This  iB  much  more  efficient  than 

computing  the  absorber  scatter  from  every  grid  point  at  each  frequency. 

A  plot  of  ceiling,  wall,  floor,  and  total  absorber  scatter  versus  frequency  is 
now  made  by  PLOTTER.  The  results  arc  also  written  to  an  output  file,  OP- 
TI0N-2.DAT.  This  file  is  shown  here: 


IT 

2.000000 

-41.61463 

-63.82076 

-10000.00 

-41.77732 

3.000000 

-46.03646 

-67.34268 

-10000.00 

-46.20884 

4.000000 

-47.63622 

-60.84136 

-10000.00 

-47.70840 

S. 000000 

-40.47342 

-61.77066 

-10000.00 

-48.73661 

6.000000 

-61.06706 

-63.36310 

-10000.00 

-61.32024 

7.000000 

-62.30607 

-64.70213 

-10000.00 

-62.65016 

8.000000 

-63.66681 

-66.86107 

-10000.00 

-63.81000 

0-000000 

-64.67887 

-66.88603 

-10000.00 

-64.84208 
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SPEW 


10.00000 

-65.40402 

-67.80018 

-10000.00 

-65.75718 

li. ooooo 

-56.32187 

-68.62801 

-10000.00 

-66.58506 

12.00000 

-67.07764 

-60.38380 

-10000.00 

-67.34083 

13.00000 

-57.77288 

-70.07003 

-10000.00 

-58.03607 

14.00000 

-58.41657 

-70.72273 

-10000.00 

-68.87076 

15.00000 

-60.01685 

-71.32108 

-10000.00 

-50.27804 

16.00000 

-60.57642 

-71.88256 

-10000.00 

-50.83861 

17.00000 

-60.10300 

-72.40814 

-10000.00 

-60.36618 

18.00000 

-60.50047 

-72.80561 

-10000.00 

-80.86268 

The  first  number  in  the  file,  17,  specifies  the  number  of  frequencies  considered. 
Each  row  then  lists  frequency  (GHz),  total  scattered  power  (dB),  scattered  power 
from  the  ceiling  (dB),  scattered  power  from  the  side  walls  (dB),  and  scattered 
power  from  the  floor  (dB).  Again,  a  value  of -10,000  dB  indicates  that  no  absorber 
scattered  power  arrived  during  the  receive  window.  CHAMBER  ends  the  run  with 
these  reminders: 

The  option  2  data  file  has  bean  arittan. 

YOU  SHOULD  PR1KT  THE  CB.L06.DiT  FILE  III 

Lastly,  let  us  select  option  3.  Again,  CHAMBER  requests  user  input.  The 
questions  asked  are  shown  (answered)  below: 

ESTER  tha  daairad  option  nunbar:  3 

Entar  tha  frequency  of  operation  (GHz):  10 

Entar  tha  taat  point  coordinata  to  be  varied  (X=1,Y=2,Z=3) :  3 

Entar  tha  starting  Z  value  (ft):  20 

Entar  tha  ending  Z  value  (ft):  28 

Entar  tha  atap  in  Z  (ft):  0.6 

Entar  tha  Teat  point  X  coordinata  (ft):  6 

Entar  tha  Taat  point  Y  coordinata  (ft):  0 


Once  the  responses  have  been  made,  CHAMBER  lists  them  and  waits  for  verifi¬ 
cation  before  continuing: 


The  frequency  of  operation  =  10.0  (GHz) 


Z  ranges  from  20.00  to  28.00  in  steps  of  0.50  (ft) 
X  =  fl.OO  (ft) 

Y  =  0.00  (ft) 

Input  OX  (t/f)  T  T 

Note  that  the  following  is  added  to  CHJLOG.DAT: 

Option  3  vas  chosen,  with  these  parameters : 


The  frequency  of  operation  =  10.0  (GHz) 

Z  ranges  from  20.00  to  28.00  in  steps  of  0.50  (ft) 

x  =  a.oo  (ft) 

Y  =  0.00  (ft) 

CHAMBER  then  computes  the  timing  data  and  all  absorber  scattering  for 
each  test  point  location.  When  working  with  each  individual  test  point,  a  message 
such  as  this  is  seen: 

The  current  test  point  is;  X  =  fl.OO  (ft) 

Y  =  0.00  (ft) 

Z  e  20.00  (ft) 

Computing  timing  data. . . 


Computing  scatter  froa  (entire)  ceiling  to  test  point... 

Coaputing  ecatter  froa  (both)  tide  ealls  to  teet  point... 

Coaputing  ecatter  froa  (entire)  floor  to  test  point... 

A  plot  of  ceiling,  wall,  floor,  and  total  absorber  scatter  versus  test  point 
location  is  now  made  by  PLOTTER.  The  results  are  also  written  to  an  output  file, 
OPT10N-3.DAT.  This  file  is  shown  here: 
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17 


20.00000 

-53.85984 

-66.60120 

-86.64442 

-54.10539 

20.50000 

-65.78187 

-66.82680 

-87.93372 

-66.14011 

21.00000 

-56.01900 

-67.16277 

-90.61380 

-56.26962 

21.60000 

-56.05049 

-67.46210 

-87.43762 

-66.37653 

22.00000 

-65.87511 

-87.64201 

-10000.00 

-66.17431 

22.50000 

-55.78619 

-67.80883 

-10000.00 

-56.06770 

23.00000 

-56.73362 

-67.92527 

-10000.00 

-58.00407 

23.50000 

-55.82059 

-67.97193 

-10000.00 

-65.88097 

* 

24.00000 

-55.49404 

-67.80015 

-10000.00 

-65.76723 

* 

24.50000 

-55.33838 

-67 ,74808 

-10000.00 

-55.59525 

26.00000 

-55.37782 

-67.80435 

-10000.00 

-55.63358 

25.50000 

-55.19563 

-67.69741 

-10000.00 

-55.44688 

26.00000 

-56.08461 

-87.58611 

-10000.00 

-55.33588 

26.50000 

-65.41229 

-67.72816 

-10000.00 

-65.87487 

27.00000 

-58.29288 

-71.97845 

-10000.00 

-68.41177 

27.60000 

-59.90308 

-10000.00 

-10000.00 

-59.90308 

28.00000 

-10000.00 

-10000.00 

-10000.00 

-10000.00 

The  first  number  in  the  file,  17,  specifies  the  number  of  test  point  locations  consid¬ 
ered.  Each  row  then  lists  the  value  (in  feet)  of  the  test  point  coordinate  that  was 
varied,  total  scattered  power  (dB),  scattered  power  from  the  ceiling  (dB),  scattered 
power  from  the  side  walls  (dB),  and  scattered  power  from  the  floor  (dB).  Again,  a 
value  of  -10,000  dB  indicates  that  no  absorber  scattered  power  arrived  during  the 
receive  window.  CHAMBER  ends  the  run  with  these  reminders: 

Tht  option  3  data  lila  hat  toaaa  Mitten. 

YOU  SHOULD  miT  THE  CELLOS. DAT  FILE  HI 


4  Before  closing  this  section,  it  is  noted  that  the  views  provided  by  the  (option  1 ) 

gray-scale  plots  of  Chapter  V  are  depicted  in  Figure  70. 

A.4  Error  Messages 

Program  CHAMBER  makes  some  simple  checks  on  its  input  data,  but  does 
not  attempt  to  detect  all  possible  errors.  However,  with  these  few  checks,  and  the 
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ample  display  of  information  provided  by  CHAMBER,  it  is  unlikely  that  bad  data 
would  go  unnoticed. 

Let  us  now  consider  error  messages  that  may  occur  while  running  program 
CHAMBER.  If  the  x1  coordinate  of  the  tips  of  the  ceiling  pyramids  is  less  than  the 
maximum  x1  value  of  the  reflector,  the  following  error  message  is  displayed: 

X  coordinato  of  tips  o 1  pyramids  on  ceiling  is  too  small ! 

If  the  y'  coordinate  of  the  tips  of  the  side  wall  pyramids  is  less  than  the  maximum 
x'  value  of  the  reflector  (which  is  also  the  maximum  y'  value  of  the  reflector),  then 
this  error  message  is  displayed: 

Y  coordinate  of  tipi  o 1  pyramids  on  aide  sail  ia  too  email! 

The  x1  coordinate  of  the  tips  of  the  floor  pyramids  should  be  less  than  zero.  If  this 
is  not  the  case,  the  following  error  message  occurs: 

X  coordinate  of  tipo  of  floor  pyramid*  ia  >  than  zero) 

In  file  CHJNPUT.DAT,  the  user  specifies  the  maximum  z1  coordinate  of  interest 
in  the  room,  and  the  z1  coordinate  of  the  target  zone  center.  Naturally,  the  target 
zone  center  should  lie  within  the  z1  extent  of  the  room.  If  not,  the  following  error 
message  occurs: 

Maximum  Z  coordinate  of  int*raat  doaa  not  go  bayond  targat 
zona  cantor  I 

CHAMBER  also  checks  that  the  test  point  is  in  the  target  region  of  the  range;  that 
is,  that  the  reflection  point  corresponding  to  the  test  point  lies  along  the  parabolic 
section  of  the  reflector.  If  this  is  not  the  case,  the  following  is  displayed; 

Test  point  not  in  targat  ragionl! 

If  any  of  these  occurs,  the  user  will  be  instructed  to  modify  the  data  in  file 
CHJNPUT.DAT,  or  to  input  different  data  from  the  terminal,  as  appropriate. 
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There  are  other  error  messages  which  will  not  occur  during  a  run  of  CHAM¬ 
BER.  However,  they  exist  in  the  routines  discussed,  and  thus  are  mentioned  for  the 
sake  of  completeness.  Recall  that  subroutine  REFL-CODE  is  a  modified  version 
of  the  reflector  code  in  [9].  This  code  computed  the  illumination  on  the  y1  >  0 
side  of  the  room  only  (since  illumination  is  symmetric  about  the  room  centerline). 
The  y 1  coordinate  of  field  points  must  be  positive.  Thus,  the  following  error  mes¬ 
sage  was  included  in  the  original  reflector  code,  and  carried  over  into  subroutine 
REFLXODE: 

Begat ive  value  of  y  la  not  allowed. 

CHAMBER  will  not  commit  this  error. 

To  understand  the  ne.it  message,  one  must  realize  that  the  original  reflector 
code  asked  the  user  to  choose  one  of  four  types  of  field  cuts.  Program  CHAMBER 
only  requires  two  types  of  cut6.  Thus,  when  subroutine  REFL-CODE  was  written, 
it  was  ensured  that  these  two  field  cuts  would  work  properly.  The  same  was  not 
done  for  the  other  two  cuts  (1  and  4).  However,  they  were  not  eliminated  from 
the  code,  in  the  chance  that  someone  may  want  to  work  with  them  in  the  future. 
For  this  reason,  if  these  cuts  are  chosen,  the  following  warning  presently  appears: 

NARVIK i  Coda  has  not  bean  nod Iliad  for  XFCUT  =1,4. 

Again,  this  will  not  occur  while  running  CHAMBER. 

Finally,  the  original  reflector  code  displays  an  error  message  and  execute 
a  FORTRAN  STOP  if  it  comes  across  a  field  point  for  which  it  cannot  find  a 
corresponding  reflection  point  (on  the  reflector).  When  subroutine  REFL-CODE 
comes  across  a  field  point  along  some  field  cut  for  which  no  reflection  point  can 
be  found,  it  will  end  its  calculations,  but  will  return  to  the  main  program  with 
its  partial  results  and  with  a  flag  set  to  indicate  what  has  happened.  The  calling 
program  should  check  this  flag  after  each  call  of  subroutine  REFL-CODE,  and 
take  whatever  action  necessary.  CHAMBER  does  check  this  flag,  and  if  set  will 
display: 

So  reflaction  point  found  -  aoaathlng  ia  arongl 
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The  message  indicates  usomething  is  wrong”  since  CHAMBER  is  specifically  de¬ 
signed  to  consider  only  field  points  that  have  corresponding  reflection  points.  This 
message  shottld  never  appear. 

One  last  point  must  be  made  on  subroutine  REFL-CODE.  While  the  mea¬ 
sured  feed  pattern  option  was  carried  over  from  the  original  reflector  code,  and 
while  CHAMBER  will  not  object  to  one  specifying  a  measured  feed  pattern  in  file 
CHJNPUT.DAT,  this  option  has  not  been  exercised  by  the  author. 
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